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COMPARATIVE PERFORMANCE ANALYSIS OF VSA AND GWO TECHNIQUES FOR 
LOAD FREQUENCY CONTROL OF SINGLE AREA FUEL CELL MICROGRID 
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Abstract. In this paper, the frequency control of a single area micro grid system is studied. 
Recently, the dynamic performance of micro-grid systems can be enhanced by effective load frequency 
control which is one of the most important aspects. The system having a fuel cell, micro turbine, diesel 
generator, battery energy storage system and flywheel energy storage system 1s modeled. To optimize 
the dynamic performance of the system, Vortex Search Algorithm (VSA) and Grey Wolf Optimizer 
(GWO) are used. Additionally, a cascade PD-PI controller is also applied to the system for enhancing 
the performance. The results are analyzed with respect to overshoot, undershoot and settling time. It can 
be clearly observed from the results that GWO algorithm provides better system performance for selected 
micro-grid system. 


Keywords: Load frequency control, micro grid, vortex search algorithm, grey wolf optimizer 
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В данной статье рассматривается частотное регулирование отдельной микросетевой 
системы. В последнее время динамические характеристики микросетевых систем могут быть 
улучшены за счет эффективного управления частоты нагрузки, что является одним из наиболее 
важных аспектов. 

Смоделирована система, имеющая топливный элемент, микротурбину, дизельный 
генератор, аккумуляторную систему накопления энергии и систему маховичного накопителя 
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энергии. Для оптимизации динамических характеристик системы используются алгоритм поиска 
вихрей (VSA) и оптимизатор серого волка (GWO). Кроме того, для повышения 
производительности в системе также применяется каскадный контроллер PD-PI. Результаты 
анализируются относительно времени перерегулирования, недорегулирования и установления. 
Из результатов ясно видно, что алгоритм GWO обеспечивает лучшую производительность 
системы для выбранной микросетевой системы. 


Ключевые слова: контроль частоты нагрузки, микросетевая система, алгоритм поиска 
вихрей, оптимизатор серого волка 


Introduction: Conventional electrical power systems are mainly consist of large number of 
components. Therefore, transmission of generated power from generation units to the over distance last 
consumer reduces the transmission efficiency [1]. In parallel, these power systems mostly use fossil fuels 
to provide as energy resources. The transfer of energy to long distances causes an increase in the fossil 
fuels used as well. 

In the globalizing world, traditional power systems are gradually being replaced by modern power 
systems [2,3]. Micro-grid systems can be considered as modern power systems. These systems usually 
contain non-fossil source (such as fuel cell, battery systems) and renewable energy sources (such as wind 
turbines, solar panels). 

As with conventional power systems, load frequency control is also a very important requirement 
for micro-grid systems. The main objective of load frequency control can be considered as follows; 


1) Keep the system frequency specific band width limits [4], 
2) Deviation of frequency made zero or closest to zero [5], 


In this study, VSA [6] and GWO [7] are applied to the selected micro-grid system to optimize load 
frequency control performances of the system. The frequency deviation obtained by both algorithms are 
compared in terms of overshoot, undershoot and settling times. Moreover, cascade PD-PD controller is 
used in these analyzes. This controller type consists of combination of PD controller and PI controller. 

This paper is organized as follows: In section 2, micro-grid model for load frequency control, VSA 
and GWO are briefly mentioned. In section 3, obtained results are illustrated. In section 4, Conclusions 
are given. 


Material and Methods 

I. Micro grid model for load frequency control: Load frequency control is essential control for 
micro grid power systems. The micro grid model used in this study is given in Figure 1. This model have 
micro turbine block, diesel generator block, battery energy storage system (BESS) block, flywheel 
energy storage system (FESS) block, micro grid (MG) system block, fuel cell, inverter block, connection 
device block and droop characteristic block. 
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Fig. 1. Micro grid model for LFC analyzes 
In this model, PD controller can be defined as follows: 
PD = Кра d Ska (1) 


Where, Кра and ka proportional and derivative gain of PD controller. PI controller can be defined 
as follows: 


ki 
PI = kpi + E (2) 
Where, kpi and ki proportional and integral gain of PI controller. Drop characteristic can be defined 
as follows: 
Drop Characteristic — - (3) 


Micro turbine can be defined as follows [8]: 


1 
1+ Тит 





(4) 


Where, Twr is time constant of micro turbine. Diesel can be defined as follows [8]: 


1 
1+sTp 





(5) 


Where, Tp is time constant of diesel. Fuel cell can be defined as follows [2]: 
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1 











1+5Тес (6) 
Where, Trc 1s time constant of fuel cell. Inverter can be defined as follows [2]: 
1 
1+51 (7) 
Where, Ti 15 time constant of inverter. Connection device can be defined as follows [2]: 
1 
1+sTcp (8) 
Where, Тсо 15 time constant of connection device. MG system can be defined as follows [2]: 
1 
(9) 





Ms+D 


Where, M is inertia constant and D is damping coefficient respectively. BESS can be defined as 
follows [8]: 


1 
1+ST BESS 





(10) 
Where, Terss 15 time constant of battery energy storage system. FESS can be defined as follows 
[2]: 


1 
1+ST FESS 





(11) 


Where, Tress 15 time constant of flywheel energy storage system. İn this study, ITSE (integral time 
square error) is used as objective function in GWO and VSA techniques. ITSE can be defined as 
following equation: 


ITSE = İ(time x (error)?)dt (12) 


II. Vortex Search Algorithm: VSO is developed by B. Dogan and T. Olmez in 2014 [6]. This 
algorithm has good balance between exploration and exploitation behavior. For this reason, the optimal 
results can be reliably found. The shape of the mixed liquids inspired the development of the VSA [6]. 

In this technique, candidate solutions are formed in the first circuit. All candidate solutions are 
controlled for the problem and the best one (iteration best) 1s determined from these and defined as global 
best solution. Founded iteration best solution 1s defined as the center of the next circuit. After then, radius 
of the circuit is decreased. New iteration best solution is compared with the previous best iteration ”s 
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solution. If new iteration best solution has minimum value rather than previous best iteration solution, 


new iteration best solution is defined as global best solution. 


This process continues until the termination criterion is met. Basically flow diagram of VSA 1s 


illustrated in Figure 2. 


Generate random solutions 


Tested all candidate solutions 


Define best iteration solution, reduce the 
radius of circuit, 








Made best iteration solution as next circuit 
center 







If best iteration solution is better than global 
best, global best= best iteration solution 


Is termination criteria met? 


YES 


Fig. 2. Flow chart of VSA 


III. Grey Wolf Optimizer: GWO is developed by S. Mirjalili et al. in 2014 [7]. This technique 
has very good capability to search of parameter space. Grey wolves are inspired the development of the 


GWO [7]. 


In this method, positions of grey wolves (candidate solutions) are formed randomly. These 


candidate solutions are search the space for obtained minimum objective function for the problem. After 


the iteration, best iteration values are determined and then the positions of candidate solution are updated. 


If the best iteration value in the next iteration 1s better than the previous best iteration value, the best 


iteration value in the next iteration 1s defined as the global best value. 


This process continues until the termination criterion is met. Basically flow diagram of GWO is 


shown in Figure 3. 
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Generate random wolves positions 


Tested all wolves according to the objective 
function 


Define best solution 


If best solution is better than global best, 
global best= best solution 


Update the wolves position 
Is termination criteria met? 
YES 










Fig. 3. Flow chart of GWO 


Results and Analysis: In this section, comparative performance results of VSA and GWO are 
discussed. Micro grid system parameters are given in Table 1. In this system, some parameters are taken 
from the reference [2], while some parameters are taken from the reference [8]. 


Table 1: Parameters of micro grid system 
R Tur Tp M D 
0.05 2 2 0.1667 0.015 


Trc Ti Tcp Твеѕѕ Tress 
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0.26 0.04 0.004 0.1 0.1 


The cascade PD-PI controller parameters with VSA and GWO are given in Table 2: 


Table 2: Controller parameters 


kpi ki Kpa Ka 
VSA -1.7945 -4.9489 4.9499 2.3019 
GWO 2.8922 4.9500 -4.7296 -1.4668 


Frequency deviations outputs are illustrated in Figure 4: 


GWO 
-----VSA 


Frequency Deviation (Hz) 





0 2 1 1.3 2 Ae 3 
Time (sec.) 
Fig. 4: Comparative frequency deviation graphic 


From this figure, it is clear that the obtained two frequency deviation outputs with VSA and GWO 
are close to each other. However, the oscillation amplitude of GWO is slightly larger compared to VSA. 
Obtained numerical overshoot, undershoot and settling time results are given in Table 3. 


Table 3: Comparative performance results 
Overshoot Undershoot Settling time (0.002 band with) 
VSA 3.411e-3 -1.258е-2 0.5054 
GWO 4.157e-3 -1.236e-2 0.4020 


Saka M., Gozde H., Eke L, Yusifbayli N., Taplamacioglu M.C.: COMPARATIVE PERFOR... 8 


It is clearly seen that from this table, although VSA gave minimum overshoot value, GWO gave 
minimum settling time and undershoot values. 


Conclusions: In this paper, load frequency control of single area having fuel cell micro grid system 
is examined. Comparative performance results performed with VSA and GWO methods. Overshoot 
value obtained by VSA approximately 9617.94 lower than GWO. However, GWO gave approximately 
261.75 lower undershoot and % 20.46 lower settling time value than VSA. It can be said that from the 
obtained results, GWO has the ability to give better performance. 
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УДК 621.311.110 


«ПРИНЦИПЫ ФОРМИРОВАНИЯ ТРЕБОВАНИЙ К ОБЪЕМАМ 
И БЫСТРОДЕЙСТВИЮ РАЗГРУЗКИ В ЭНЕРГОРАЙОНАХ ПРИ 
ДЕЙСТВИИ ДЕЛИТЕЛЬНОЙ АВТОМАТИКИ» 
«PRINCIPLES BEHIND THE REQUIREMENTS ТО THE EXTENT AND 
SPEED OF LOAD-SHEDDING BY ISLANDING AUTOMATIONS IN 
MICROGRIDS» 


Илюшин Павел Владимирович!, Куликов Александр Леонидович” 
Ilyushin Pavel Vladimirovich, Kulikov Aleksandr Leonidovich 
"Институт энергетических исследований Российской академии наук, 
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Energy Research Institute of the Russian Academy of Sciences, Moscow, Russia 
“Нижегородский государственный технический университет 
им. Р.Е. Алексеева, г. Нижний Новгород, Россия 
Nizhny Novgorod State Technical University n.a. К.Е. Alekseev, 

Nizhny Novgorod, Russia 


Аннотация: Анализ статистических данных показывает, что выделение частей 
энергосистем (энергорайонов) в основной режим происходит, как правило, в результате 
аварийных возмущений в ремонтных схемах сети. Применение делительной автоматики 
позволяет осуществить выделение электростанции (генерирующей установки) или 
энергорайона с целью сохранения в работе потребителей собственных нужд, а также 
максимального числа электроприемников сторонних потребителей, подключенных к 
сетям среднего напряжения электростанции. Показано, что выделение энергорайона 
может сопровождаться как понижением частоты, если устойчивость электродвигателей в 
нагрузке не нарушается, так и ее повышением, если возникающий дефицит мощности 
провоцирует возникновение лавины напряжения. Результаты расчетов переходных 
процессов показывают, что величина критического дефицита мощности существенно 
зависит от наличия инициирующего выделение КЗ и состава нагрузки, а также 
распределения реактивной нагрузки в доаварийном режиме между электростанцией 
энергорайона и энергосистемой. Обосновано, что устойчивость нагрузки является 
определяющей при выборе объемов и быстродействия разгрузки в выделяемом 
энергорайоне, что требует корректного учета состава и параметров нагрузки при 
моделировании. Показано, что быстродействие разгрузки особенно критично, если 
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электродвигатели в результате близкого КЗ быстро тормозятся, что может спровоцировать 
возникновение лавины напряжения и нарушение электроснабжения всех потребителей. 


Abstract: Statistics shows that microgrids are islanded when an emergency-related 
disturbance occurs whilst under maintenance. Islanding automations help isolate a power plant 
(generator set) or a microgrid to keep online the inhouse loads and as many third-party consumers 
connected to the medium-voltage grid of the plant. Apparently, islanding a microgrid may cause 
the frequency to drop if the motors constituting the load are not destabilized; or to rise, if the power 
deficit triggers a voltage collapse. Transient calculations show that the critical power deficit value 
depends on whether there is an islanding-initiating short-circuit, on the load composition, as well 
as on the pre-emergency distribution of reactive load between the microgrid's power plant and the 
general power system. The paper shows that the load stability is the decisive factor when it comes 
to configuring the extent and speed of load-shedding in an islanded microgrid; this requires due 
account of the composition and parameters of the load in modeling. The paper also shows that the 
speed of load-shedding is most critical if a nearby fault causes motors to decelerate quickly, which 
might trigger a voltage collapse and thus a power outage for all consumers. 

Ключевые слова: делительная автоматика, электростанция, энергорайон, островной 
режим, дефицит мощности, нарушение устойчивости электродвигателей, лавина 
напряжения, быстродействие разгрузки. 

Keywords: islanding automation, power plant, microgrid, islanding, power deficit, motor 
destabilization, voltage collapse, load-shedding speed. 


Введение 


В последние годы в условиях цифровизации электроэнергетических систем наблюдается 
тенденция к росту величины износа генерирующего и электросетевого оборудования в 
Единой энергетической системе России (ЕЭС России), которая наиболее ярко проявляется 
в распределительных сетях. Это обусловлено как объективными причинами (сдерживание 
роста тарифа на услуги по передаче электрической энергии по сетям; снижение платежной 
дисциплины; рост величины не технологических потерь и др.), так и субъективными 
(большие потери электроэнергии в сетях — до 20%; высокие операционные издержки B 
субъектах электроэнергетики; низкая эффективность капитальных вложений и др.) [1-3]. 


В этих условиях растет количество аварий, в том числе системных, с тяжелыми 
технологическими последствиями и значительными ущербами у промышленных 
потребителей [4-6]. 


В результате возникновения аварий, в том числе в ремонтных схемах сети, с достаточной 
регулярностью происходят выделения частей энергосистем (энергорайонов), находящихся 
в зоне централизованного электроснабжения, в основной режим работы. 


Под энергорайоном понимается совокупность объектов энергосистемы, расположенных 


Ilyushin Pavel V., Kulikov Aleksandr L.: PRINCIPLES BEHIND THE... 12 


на части обслуживаемой ею территории. 


Под островным режимом подразумевается режим работы энергорайона с одной или 
несколькими электростанциями, каждая из которых может иметь в своем составе одну или 
несколько генерирующих установок (ГУ), и нагрузку, который допустим по всем 
условиям электроснабжения и электропотребления, образующийся при отключении линий 
электропередачи (ЛЭП), связывающих энергорайон с энергосистемой (в результате КЗ; 
без КЗ) и существующий до момента его синхронизации с энергосистемой. 


Учитывая технологические особенности применяемых на электростанциях ГУ и 
технологических линий на промышленных предприятиях обеспечить их надежное 
функционирование без выделения в островной режим, в случае возникновения (развития) 
аварии в энергосистеме, в большинстве схемно-режимных ситуаций, невозможно. Для 
этих целей на электростанциях (в энергорайонах) применяется делительная автоматика 


(ДА) [7]. 


Под делительной автоматикой понимается особый вид противоаварийной автоматики 
(ПА), осуществляющей выделение в аварийных схемно-режимных условиях 
электростанции (генерирующей установки) или энергорайона от энергосистемы с целью 
сохранения в работе потребителей собственных нужд, а также максимального числа 
электроприемников сторонних потребителей, подключенных к сетям среднего 
напряжения электростанции [8]. 


Делительная автоматика, обеспечивающая выделение электростанции (отдельных 
генераторов) для питания потребителей собственных нужд при возникновении аварийного 
режима в энергосистеме известна, в основном, в виде частотной делительной автоматики 
(ЧДА). Данная автоматика срабатывает при снижении частоты в энергосистеме до 46—47,5 
Гц. Ее задача состоит в том, чтобы при возникновении в какой-либо отделившейся от ЕЭС 
России части со значительным дефицитом мощности при недостаточном объеме нагрузки, 
заведенном под действие устройств автоматической частотной разгрузки (АЧР), 
предотвратить полное погашение электростанции из-за недопустимого снижения частоты, 
сохранив электроснабжение потребителей собственных нужд [9]. 


Целью статьи является анализ особенностей переходных процессов при выделении 
электростанции (энергорайона) в островной режим без КЗ и в результате КЗ, оценка 
влияния особенностей прилегающей сети и параметров нагрузки на характеристики 
переходных процессов, а также разработка принципов формирования требований к 
объемам и быстродействию разгрузки при действии ДА. 
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Результаты расчетов переходных процессов при выделении знергорайона c Ш ОС в 
островной режим 

На рис. 1 представлена упрощенная однолинейная схема парогазовой электростанции 
(ПГЭС), содержащей четыре газотурбинные установки (ГТУ) в одновальном исполнении 
(ГТ-Ти ГТ-2 — по 43 МВт, ГТ-4 и ГТ-5 — по 45 МВт) и две паротурбинных установки 
(ПТУ) — (IIT-3 и ПТ-6 по 30 МВт) в составе двух парогазовых установок (ГТ-1 + ГТ-2 — 
ПТ-3 и ГТ-4 + ГТ-5 — ПТ-6). 


кПС-1 110 кВ k [1C-2 
ЛЭП-4 | ЛЭП-3 ; nən-1 ЛЭП-2 
1 m 





PYCH-10 KB PYCH-10 KB 


Рис. 1. Упрощенная однолинейная схема Ш ЗС. 


Рассмотрим энергорайон с суммарной генерацией 236 МВт (при всех включенных ГУ на 
Ш ЗС) и суммарной нагрузкой 279 МВт. В данном энергорайоне проектом планируется 
установка ДА. Приведенные в статье результаты расчетов переходных процессов 
выполнены в реальной схеме сети в ПК Мустанг. 


Проведем анализ переходных процессов в энергорайоне c Ш ЗС при различных 
дефицитах мощности в случае отключения ЛЭП-1 и ЛЭП-2 или одной из них, 
связывающих энергорайон с ПС-2 при выведенной в ремонт второй. На момент 
выполнения расчетов режимов ЛЭП-3 и ЛЭП-4, связывающие энергорайон с ПС-1, 
находились в стадии строительства и не были введены в работу. В расчетах действие 
устройств ПА в энергорайоне не учитывалось с целью выявления влияния параметров 
нагрузки и предаварийного режима на характеристики переходных процессов. 


В серии расчетов учитывались различные величины дефицита мощности в энергорайоне, 
который рассчитывался по формуле: 
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р=|1- 2^ | 100% 


2.Р, 


где: У Р. – суммарная располагаемая величина генерации B знергорайоне; > Py – суммарная 
величина нагрузки в энергорайоне в предаварийном режиме. 


При выполнении расчетов варьировалось значение дефицита мощности. Было принято 
допущение, что величина нагрузки в сети 110 кВ не изменялась. Рассматривалось 
увеличение нагрузки в КРУ-20 кВ Ш ЗС в диапазоне от 0 до 200 МВт. Если критическое 
значение дефицита мощности (Dkp) оказывалось меньше соответствующего расчетной 
нагрузке в энергорайоне, то снижались все нагрузки на один и тот же процент. В 
рассматриваемой схеме превышение значения Dkp приводит к тому, что переходный 
процесс со снижением частоты заканчивается лавиной напряжения и сопровождается 
значительным сбросом нагрузки с последующим повышением частоты [10, 11]. 


С учетом исходных данных в рассматриваемом энергорайоне Dy, = 43,3%. На рис. 2a-e 
приведены результаты расчетов переходных процессов при выделении энергорайона с 
дефицитом мощности меньшим и большим величины Ок. Основные технические 
характеристики эквивалентных АД приведены в таблице 1. 


Таблица 1. Технические характеристики эквивалентных АД. 


Коэффициент 
загрузки АД 


Собственные нужды ПГЭС 


Промышленная нагрузка на 
напряжении 10 и 20 кВ 


Эквивалентная нагрузка в сети 110 
кВ 


(зимний максимум) 


Эквивалентная нагрузка в сети 110 
кВ 


(летний максимум) 





EEEC Volume 11 Nel 2021 


Частота - f (Гу) Напряжение - Ш/ ном (0.e.) 


54 | 12 
10 
| 0,9 
0,8 
0,7 
0,6 
0,5 
0,4 
0,3 


0,2 





Относительные скольжения AJ] - s (96) 
100 


АД на стороне 10 кВ ПС-2 


Остальные АД к энергорайоне 


Lc 


15 
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Активные мощности ГУ FTT 3C - Р (МВт) 


80 
72 


64 











6) 


Рис. 2. Переходные процессы в сети 110 кВ энергорайона после срабатывания ДА при 
D < Dep = 43,3% и D > Diş: a – частота и напряжение в сети 110 KB; 6 – относительные 
скольжения АД; в – активные мощности ГУ ПГЭС. 


Анализ рис. 26 позволяет сделать вывод, что лавину напряжения в энергорайоне 
провоцирует нарушение устойчивости части асинхронных двигателей (АД). Наименьший 
запас устойчивости зафиксирован у АД, подключенных к шинам 10 кВ ПС-2. 

На рис. 26 часть АД при D = 45% полностью останавливаются и отключаются 
действием устройств релейной защиты (РЗ), а у части начинаются процессы самозапуска. 
Самозапуски АД могли бы закончиться успешно, если бы напряжение в сети 110 кВ не 
снизилось ниже 60% or (һом. В реальной схемно-режимной ситуации самозапуск будет 
неуспешным и АД будут отключены действием устройств P3 по снижению напряжения или 
технологическими защитами приводимых механизмов [12-14]. 

Анализ приведенных расчетов переходных процессов позволяет сделать вывод, что при 
одном и том же исходном режиме изменения параметров нагрузки приводят к изменению 
величины Dkp, как показано в таблице 2. 
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Таблица 2. Изменение О в зависимости OT параметров нагрузки. 


Изменяемые параметры нагрузки 
Исходное состояние расчетной схемы 


Реактивные нагрузки увеличены в 1,5 раза со средним созф = 0,831 (вместо 


0,913) 
Доля АД в непромышленной нагрузке увеличена до 60% (вместо 40%) 


Отключена одна паротурбинная ГУ (ПТ-3) 


При изменении коэффициента мощности нагрузки (COS@y) изменение Dy, не существенно, 
а в случае увеличения доли АД в составе нагрузки Ha 20% величина Dkp снижается почти в 
1,6 раза, что требует при выполнении расчетов уделять особое внимание корректности 
моделирования нагрузки в сети. Еще большую роль в снижении величины Р играет 





снижение величины располагаемой мощности в рассматриваемом энергорайоне. В случае 
отключения одной паротурбинной ГУ на ПГЭС [хр снижается почти в 2 раза. 


При выполнении расчетов был выявлен еще один фактор, оказывающий существенное 
влияние на величину Dkp, как показано на рис. 3. Этим фактором является напряжение 
( (о) в доаварийном режиме (перед выделением) на шинах 110 кВ ПС-2, которая 
определяется режимом работы ПС 500/220/110 кВ. 


40 


30 


20 





Критический дефицит мощности - О (90 


10 
104 108 112 116 120 


Напряжение на шинах 110 кВ ПС-2 
в доаварийном режиме - Uo (кВ) 


Рис. 3. График зависимости Dp от Uo. 


Используемые в расчетах величины Uo относятся к доаварийным режимам работы 
прилегающей сети и характеризуют ее особенности, при этом нагрузка в энергорайоне 
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соответствует расчетной схеме, a D = 24,6%. 


Анализ рис. 3 позволяет выявить скачкообразное изменение Dep при незначительных 
изменениях величины Uo в диапазоне от 112 до 116 кВ. Данный эффект обусловлен 
изменением загрузки ГУ ПГЭС по реактивной мощности, как показано на рис. 4. Это 
происходить в случае, если на смежной ПС 500/220/110 кВ повышается напряжение на 
шинах 220 кВ. При этом ГУ на ПГЭС в рассматриваемом энергорайоне разгружаются по 
реактивной мощности в пределах допустимого диапазона (0-133,5 МВАр). 


(МВАр) 
> 


~ 
ж- 





104 108 112 116 120 
Напряжение на шинах 110 кВ ПС-2 B 
доаварийном режиме - Uo (кВ) 


Суммарная выработка реактивной мощности 


ПГӘС в доаварийном режиме - Опгэ‹ 
I 
m 


Puc. 4. График изменения выработки Опгэс в зависимости от величины Uo 
в доаварийном режиме. 


Таким образом, чем меньшую реактивную мощность вырабатывают ГУ ПГӘС (Опгэс) в 
доаварийном режиме, тем глубже проваливается напряжение в начальный момент 
времени при выделении энергорайона в островной режим. Кроме того, в этом случае 
увеличивается общее время переходного процесса, что обусловлено временем, 
необходимым для реализации форсировки возбуждения ГУ Ш ЗС (нарастания токов 
возбуждения в генераторах до потолочных значений). В рассматриваемом режиме 
снижается запас устойчивости АД и существенно понижается значение Ок, 115-171. 


В режиме, когда Опгәс = 0 (например, при Uo > 117 кВ) незначительное повышение Uo 
приводит к повышению токов возбуждения ГУ, пропорциональному значениям Uo для 
обеспечения одного и то же значения Опгэс = Qnin. При этом снижаются скольжения АД в 
доаварийном режиме и несколько повышается Dy». Переходные процессы при четырех 
разных значениях Uo и D = 20,3% (точки 1, 2, 3, 4 на рис. 3) показаны Ha рис. 50-0. 
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Активные мощности ГУ ПГЭС - P (MBr) 


60 


54 





Реактивные мощности ГУ ПГЭС - Q (MBAp) 
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Относительные скольжения AJ] - s (96) 
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Рис. 5. Переходные процессы после срабатывания ДА при четырех разных значениях Uo и 
D = 20,3%: a — частота; б – напряжение в сети 110 кВ; e — активные мощности ГУ Ш ЗС; 
г — реактивные мощности ГУ Ш ЗС; д — относительные скольжения АД. 


Рассмотрим подробнее левую часть графика зависимости Dkp (Uo) на рис. 3. Режимы с 
большими значениями D (при Uo < 112 кВ) приводят к увеличению потерь напряжения B 
сети и снижению запасов устойчивости АД. В этих условиях область допустимых 
дефицитов мощности (D < Dkp) ограничена устойчивостью АД, которая дополнительно 
понижается при снижении Uo. 


Пример переходного процесса с нарушением условий существования режима при D = 50% 
> Dkp в условиях низкой величины Uo (точка 5 на рис. 3) приведен на рис. 6 
(устанавливается f = 51,65 Гц; U < 45 кВ в сети 110 KB). 


Следовательно, величины напряжения на шинах Ш ЗС в рассмотренных схемно- 
режимных условиях должны находиться в существенно более узком коридоре, чем 
регламентированные значения в нормативно-технических документах для нормальных (+ 
5%) и послеаварийных (+ 10%) режимов. 
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Относительные скольжения AJ] - s (96) 
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Рис. 6. Переходный процесс при действии ДА при Uo= 109,28 KB и D = 50%. 


Определение уровней напряжения в контрольных точках сети и их поддержание при 
ведении режимов должно основываться на проведении многовариантных расчетов 
переходных процессов для предотвращения нарушений устойчивости нагрузки и 
возникновения лавины напряжения. 


Особенности переходных процессов при выделении энергорайона c ПГЭС в 
результате КЗ 
Рассмотрим схемно-режимную ситуацию, когда ЛЭП-1, отходящая к ПС-2 отключена, а 
на ЛЭП-2 возникает трехфазное КЗ вблизи шин ПС-2. 


Результаты расчетов переходных процессов существенно отличаются от рассмотренных 
выше, где отключение ЛЭП, связывающей ПС-2 с ПС 500/220/110 кВ приводило к 
понижению частоты в результате дефицита мощности, при уровнях напряжения в узлах, 
близких к нормальным значениям. 


Рассмотрим режим, в котором ГУ ПГЭС работают с выдачей Оном, а при выделении 
знергорайона в результате трехфазного КЗ с D = 24,6% (puc. 3) устойчивость АД в 
нагрузке нарушается даже при {кз = 0,06 с (puc. 7а-г). 
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Рис. 7. Переходные процессы после выделения энергорайона в результате КЗ с D = 24,6% 
при /кз = 0,06 и /кз = 0,04 с: a – частота и напряжение в сети 110 кВ; 6 – относительные 
скольжения АД; в – активные мощности ГУ ПГЭС; г – реактивные мощности ГУ ПГЭС. 
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Анализ рис. Та-г позволяет выявить два принципиальных аспекта. 


- если отключается ЛЭП, связывающая ПС-2 с ПС 500/220/110 кВ 
(энергосистемой), то единственной электростанцией в энергорайоне будет ПГЭС, 
расположенная на его краю с суммарной генерируемой мощностью меньшей, чем 
суммарная нагрузка. В этом случае запас устойчивости электродвигателей в нагрузке 
низкий и при КЗ велика вероятность возникновения в энергорайоне лавины напряжения; 

— вероятность возникновения B энергорайоне лавины напряжения в основном 
зависит от состава и параметров нагрузки. В рассматриваемом энергорайоне значительная 
часть мощности потребляется системами кондиционирования и вентиляции высотных 
административных зданий (электродвигатели имеют момент сопротивления, близкий к 
независимому от скорости вращения — низкий запас устойчивости). 

Анализ переходных процессов в энергорайоне c Ш ЗС показал, что имеются два участка 
сети, где запасы устойчивости электродвигателей в нагрузке наименьшие. Поэтому, если 
из-за внешних факторов возникает значительное снижение напряжения, то 
опрокидывания АД начинаются именно на этих участках сети. При этом в остальных 
частях энергорайона либо И > кр, либо возникает лавина напряжения, которая 
охватывает весь энергорайон. Следовательно, разгрузку сети при действии ДА следует 
выполнять на тех участках, где запасы устойчивости электродвигателей в нагрузке 
наименьшие. 


Конкретные выводы о возможности возникновения в том или ином энергорайоне лавины 
напряжения с нарушением работы значительной части электроприемников можно делать 
только на основе детального изучения параметров нагрузки и переходных процессов с ее 
участием. 


Важно отметить, что выделение энергорайона в результате КЗ существенно повышает 
вероятность глубокого снижения напряжения без снижения частоты с нарушением работы 
значительной части электроприемников. Данный фактор следует учитывать при выборе 
алгоритмов работы и параметров настройки ДА и других устройств ПА в энергорайоне. 


Выбор мест реализации управляющих воздействий от ДА или других устройств ПА, 
осуществляющих разгрузку, зависит от состава и параметров нагрузки. Их правильный 
выбор позволяет предотвратить нарушение работы всего энергорайона при возникновении 
дефицита мощности, что требует уделять особое внимание корректности моделирования 
нагрузок [18]. 


Принципы формирования требований к быстродействию разгрузки 


Вопрос о быстродействии разгрузки (реализации управляющих воздействий от ДА на 
отключение нагрузки — ОН) особенно актуален не в тех схемно-режимных условиях, когда 
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U > О, а частота снижается, но B тех случаях, когда возможно быстрое развитие лавины 
напряжения [19]. 


В настоящее время реально обеспечить быстродействие разгрузки (реализации ОН), 
которое соответствует запаздыванию от момента срабатывания пускового органа ДА до 
момента отключения электроприемников, At = 0,6 c [20]. В условиях, когда выделение 
энергорайона не сопровождается КЗ возможно обеспечить стабилизацию параметров 
режима в энергорайоне, так как скольжения АД малы. Фактическому отключению 
нагрузки соответствует момент времени: 


Гон = Аѓпо + At + їкз, 


где Atno — время срабатывания пускового органа ДА; /кз = 0,2 c — длительность КЗ, 
ликвидируемого устройствами P3. 


При возникновении КЗ электродвигатели в энергорайоне начинают быстро тормозиться и 
к моменту разгрузки их максимальные скольжения могут превышать 60-70%. В этом 
случае без реализации ОН в больших объемах стабилизировать параметры режима не 
представляется возможным [21 ]. 


Рассмотрим аварийный режим, возникающий в результате трехфазного КЗ на ЛЭП-1 
рядом c ПС-2 c щз= 0,2 с при различном быстродействии разгрузки: ѓон = 0,2 с, юн = 0,6 с 
и Гон = 1,0 с. Результаты расчетов переходных процессов представлены на рис. 8а-в 
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Рис. 8. Переходные процессы при выделении энергорайона с ПГЭС в результате 
трехфазного КЗ с їкз= 0,2 с при объеме OH 33,5% в различным быстродействием (fon = 0,2 
с, ton = 0,6 си ѓон = 1,0 с): a – частота и напряжение в сети 110 KB; 6 – активные мощности 

ГУПГЭС; в — относительные скольжения АД. 


Анализ рис. 8 позволяет сделать вывод, что при больших объемах ОН ее быстродействие 
не играет решающей роли, однако это приводит к нарушению электроснабжения 
значительного количества электроприемников. Поэтому, быстродействующая разгрузка 
должна применяться при проектировании ДА с целью минимизации объемов ОН в случае 
угрозы возникновения лавины напряжения в островном режиме работы энергорайона 
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Выводы 


Выделение электростанции (энергорайона) действием делительной автоматики с 
потребителями собственных нужд и электроприемниками сторонних потребителей в 
случаях, когда величина нагрузки в энергорайоне превышает располагаемую мощность 
питающих ее генерирующих установок, может сопровождаться понижением частоты, если 
устойчивость электродвигателей в нагрузке не нарушается, или ее повышением, если 
возникающий дефицит мощности провоцирует лавину напряжения. 


При выборе управляющих воздействий делительной автоматики, которые следует 
реализовывать при выделении электростанции (энергорайона), требуется определить, 
возможны ли схемно-режимные ситуации, при которых нарушения устойчивости 
электродвигателей приводят к возникновению лавины напряжения в энергорайоне. Если 
это реально, то необходимо предусматривать управляющие воздействия, отключающие 
необходимый объем нагрузки по факту выделения энергорайона, обеспечивая 
поддержание напряжения на шинах нагрузки выше критических значений, с целью 
сохранения условий электроснабжения электроприемников, которые остались в работе. 
В случае выделения электростанции (энергорайона) в результате КЗ значительно 
повышается вероятность того, что это выделение приведет к глубокому снижению 
напряжения, без снижения частоты, и последующему нарушению работы значительной 
части электродвигателей в нагрузке. 

Величина критического дефицита мощности зависит от наличия инициирующего 
выделение КЗ и состава нагрузки, но в большей степени от распределения реактивной 
нагрузки в доаварийном режиме между электростанцией энергорайона и энергосистемой. 
Расчетами обосновано, что чем меньше выработка реактивной мощности генерирующими 
установками в энергорайоне, тем больше вероятность того, что дефицит мощности 
окажется больше критического и спровоцирует лавину напряжения. 

Объемы разгрузки при действии делительной автоматики должны определяться на 
основании результатов многовариантных расчетов переходных процессов, что требует 
корректного моделирования нагрузки в выделяемом энергорайоне. От доли 
электродвигателей в составе суммарной нагрузки решающим образом зависит уровень 
устойчивости нагрузки в целом. 

Быстродействие разгрузки особенно критично в тех схемно-режимных условиях, когда 
электродвигатели начинают быстро тормозиться (в результате близкого КЗ) и к моменту 
реализации разгрузки их скольжения превышают 60-70%, что может привести к быстрому 
развитию лавины напряжения. 
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Assessment of peneatration level of variable renewable generations 
on the flexibility of power system 


Yusifbayli Nurali!, Nasibov Valeh”, Guliyev Huseynqulu'Aghaliyev Nicat! 


! Azerbaijan Technical University, Baku, Azerbaijan 
? Azerbaijan Scientific Research and Design-Research Energy Institute, Baku, Azerbaijan 


Annotation: The study evaluates the impact of the level of penetration of variubale 
renewables generation sources on the flexibility of power systems on different phases in the context 
of the transformation of power systems to green technologies. The impact of solar and wind turbines, 
which make up the bulk of renewable energy sources, on the operating modes of power systems was 
also considered, as well as the forecast values of the equivalent inertia constant of the system to 
prepare the Azerbaijani power system for such a transformation. 

Key words: Varuable renewbales sourses, flexibility, penetration faze, power system, inertia 
constant. 


Doyison borpa olunan generasiya monbolorinin inteqrasiya soviyyosinin 
EES-nin cevikliyino tosirinin qiymotlondirilmosi 
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Annotasiya: İşdə elektroenegetika sistemlərinin yaşıl texnonogiyalara transformasiyası 
şəraitində dəyişən bərpa olunan əsaslı generasiya mənbələrinin penetrasiya səviyyəsinin müxtəlif 
fazaları üzrə enerji sistemlərinin çevikliyinə təsiri qiymətləndirilmişdir. Həmçinin, dəyişən bərpa 
olunan enerji mənbələrinin əsas hissəsini təşkil edən günəş və külək qurğularının enerji sistemlərin iş 
rejimlorino təsiri ayrılıqda baxılmış, eləcə də Azərbaycan elektroenergetika sisteminin belə 
transformasiyaya hazır olması üçün sistemin ekvivalent inersiya sabitinin proqnoz qiymətləri 
araşdırılmış və kompleks təkliflər hazırlanmışdır. 

Açar sözlər: Dəyişən bərpa olunan mənbələr, çeviklik, nüfuz fazası, güc sistemi, ətalət 
sabitidir. 


Giriş 

Elektroenergetika sektorunun yenidən qurulması üçün geniş imkanlar açan bərpa olunan enerji 
mənşəli (əsaslı) generasiyanın “capex” və “opex” üzrə azalan xərcləri, rəqəmsal texnologiyaların 
inkişafı və bu mənbələr hesabına istehsal olunan enerji enerji təhlükəsizliyinin təmin olunmasında 
mühüm rol oynamaqdadır [1]. Külək (KQ) və günəş fotovoltaik (FVQ) qurğuları 2040-cı ilo qədər 
əlavə elektrik enerjisi istehsalının yarıdan çoxunu və BMT-nin İDÇK (TK 21, Paris Sazişi) üzrə iqlim 
doyisikliyi öhdəliklərini yerinə yetirmək baxımından demək olar ki, bütün artımı təmin edəcəkdir. 
Hökumətlər və tənzimləyicilər texnoloji dəyişiklik tempi vo enerji sistemlərinin çevik istismarına 
artan ehtiyacla ayaqlaşmaq üçün sürətli hərəkət etməli olacaqlar. Eneryi akumlyasiya (saxlanc) 
texnologiyaları bazarının surətli inkişafı, elektrik nəqliyyat vasitələri və şəbəkə arasında interfeys, 
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məlumatların məxfiliyi və s. kimi məsələlər istehlakçılar və milli elektrik enerjisi sistemləri üçün 
potensial risk mənbəyinə çevrilməkdədir [2]. 

Son onilliklərin dinamikasına görə elektrik enerjisi istehsalında artım ümumi enerji 
tələbatından daha sürətlidir. Beynəlxalq Energetika Agentliyinin (BEA) hesablamaları göstərir ki, 
elektrik enerjisindon istifadə ümumi enerjiyo olan tələbatı iki dəfədən çox üstələyir ki, bu da onun 
müasir dövrdə ölkələrin iqtisadiyyatında əhəmiyyətli rolunu bir daha təsdiqləyir. Bu artımın 
arxasında çox saylı amillər dayanır — bunlardan elektrik avtomobilləri və bərpa olunan enerji 
mənbələrindən (BEM) birbaşa istifadə özünü xüsusilə qabarıq büruzə verir. Hazırda, enerjinin son 
istehlakında neftin yarısından az olan elektrik enerjisinin payı 2040-ci ilə qədər qlobal səviyyədə nefti 
üstələyə bilər. Exxon Mobil şirkətinin proqnozlarına görə, qlobal miqyasda elektrik eneryisinə olan 
tələbatın 6076 artacağı gözlənilir. FV sistemləri qlobal miqyasda quraşdırılmış elektrik istehsalı 
gücünün ən böyük komponenti olacaqıdr. 2040-cı По qədər aşağı karbonlu mənbələrin ümumi elektrik 
enerjisi istehsalının yarıdan çoxunu təmin edəcəyi gözlənilir. Külək və günəş fotovoltaik sistemləri 
aparıcı mənbələrə çevriləcəkdir, lakin hidroenerji (2040-cı ildə ümumi istehsalın 1596-1) və nüvə (8%) 
enerjisi də óz əsas paylarını saxlayacaqdır. Bərpa Olunan Enerj Mənbələri üzrə Beynəlxalq 
Agentliyin (BEMBA) "2050-cı ilə qədər Yol Xəritəsi" sənədinə əsasən elektrik təchizatı yalnız BEM- 
dən ucuz qiymətə alına bilən “təmiz” elektrik enerjisi ilə birlikdə emissiyaların azaldılması üçün əsas 
həll yolu kimi ortaya çıxa bilər. Ümumi enerii istehlakında elektrik enerjisinin payı qlobal səviyyədə 
2050-cı ilə qədər 50%-э qədər (hazırda 20%) artmalıdır. Bu halda BEM enerji istehlakının üçdə 
ikisini və elektrik enerjisi istehsalının 86%-ni təşkil edəcəkdir. Kütləvi elektriklosdirmo ilo iki 
dəfədən artıq artan bərpa olunan eneryi əsaslı elektrik enerjisi CO» emissiyalarını 60% azaltmaqla 
energetika sektorunda zəruri azalma miqdarında ən böyük paya sahib olacaqdır. Belə tendensiya, 
xüsusilə bərpa olunan enerji mənşəli (əsaslı) elektrik eneryisinin daha rəqabətli olacağı və sabaka 
infrastrukturunun bu gücləri qəbul etməyə hazır olacağı 2025-ci ildən sonrakı dövr üçün nəzərə 
alınması vacibdir. 

Bundan əlavə, iqlim dəyişikliyi və ətraf mühit, o cümlədən Paris Sazisi çərçivəsində dövlətlərin 
öhdəlikləri ilə əlaqəli artan tədbirlər və qlobal miqyasda karbohidrogen yanacaq aktivlərindən 
uzaqlaşma tendensiyası bütün Dünyada BEM-ə və effektiv enerjiyo transformasiyanin daha da 
sürətlənməsini gözləməyə əsas verir. BEA, enerji sektorunda bu transformasiyanın həm də bir fürsət 
olmasını vurğulayır. Bu transformasiya daha sürətli iqtisadi böyüməyə, daha çox is yerinin açılmasına 
və ümumi sosial rifahın yaxşılaşmasına zəmin yarada və nəticədə, 2050-ci ilə qədər transformasiyanın 
illik əlavə xərclərindən 3-7 dəfə çox olan vəsaitə qənaət etmək mümkün ola bilər. 2050-cı ila qədər 
епегјі transformasiyası cari vəziyyətin saxlanması ssenarisi ilə müqayisədə ÜDM-in 2,5% 
çoxalmasını və qlobal məşğulluğun 0,2% artmasını təmin edə bilər. Adekvat və ədalətli 
transformasiyani təmin edəcək siyasətin qəbul olunması müxtəlif ölkələr, bölgələr və icmalar üçün 
faydaları artıracaq və sosial bərabərsizlikləri aradan qaldıracaqdır. Qlobal energetika sistemini 
dəyişdirməklə elektrik enerjisi bütün Dünya əhalisi üçün əlçatan olacaq və enerji təhlükəsizliyi 
yaxşılaşacaqdır. 

Enerji sektorundakı qlobal vəziyyət, qlobal istixana qazlarl (İQ) emissiyalarında ciddi 
azalmalara ehtiyac olduğunu tövsiyə edən son elmi dəlillərlə xarakterizə olunsa da, digər tərəfdən 
statistik məlumatlar enerji ilo əlaqəli emissiyaların 2018-ci ildə tarixi pik həddə çatdığını göstərir. 
BEM-lərin lokomotiv olduğu enerji modellərinə surətli transformasiya gözləntiləri milli 
strategiyalara bir anda təsir etmir, belə ki, dövrümüzün karbohidrogen yanacaqlardan çox asılı olan 
enerji sistemlərində mövcud olan reallığı köklü surətdə dəyişmiş olur. BEA-nin (Noyabr 2019) dərc 
etdiyi “Dünya Enerji İcmalı” bu vəziyyəti bütün inkişaf etməkdə olan və inkişaf etmiş ölkələrdə 
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müxtəlif mümkün perspektivləri (2040-a qədər) vo enerji siyasəti qərarlarında tələb olunan 
dəyişiklikləri araşdıran fərqli ssenarilər kontekstində təqdim edir. Bundan əlavə, BEMBA-nın 
hesablamalarına görə qlobal İQ emissiyaları indidən başlayaraq 2050-ci ilə qədər ildə 3,596 
azaldılmalı və bu tendensiya həmin müddətdən sonra da davam etməlidir. 

Yuxarıda qeyd olunan global tendensiyalar inkişaf etməkdə olan Azərbaycan Respublikasına 
da samil oluna bilər. Fridrix Ebert Fondunun "Cənub-Şərqi və Şərqi Avropada, Cənubi Qafqazda vo 
Orta Asıyada Enerji Transformasiyası" adlı tədqiqatına əsasən Cənubi Qafqaz regionunda bərpa 
olunan enerji texnologiyalarının tətbiqi vo enerji səmərəliyi tədbirlərinin icrası çox yüksək texniki 
potensiala malik olmasına baxmayaraq olduqca yavaş sürətlə gedir. Azərbaycanla bağlı son bir neçə 
tədqiqat bu açıqlamanı təsdiqləyir və göstərir ki, ölkədə enerji sistemlərinin karbonsuzlasdirilmasi 
texniki cəhətdən mümkündür və sosial, iqtisadi və ekoloyi baxımdan faydalı hesab olunur. 


1. BEM-in mövcud vəziyyəti və proqnozlar 


Dəyişən bərpa olunan enerji mənbələri (DBEM) əsaslı generasiya mənbələrinin geniş tətbiqi 
fonunda elektroenergetika sistemlərinin (EES) çevikliyinin və enerfi təhlükəsizliyinin əhəmiyyətini 
kifayət qədər artırmışdır. Eyni zamanda, kiber təhlükəsizlikdən tutmuş kəskin dəyişkən hava şəraitinə 
qədər yeni təhlükələr hökumətlərdən davamlı diqqət tələb edir. Aparılan hesablamalara görə 2018-ci 
ildə qlobal eneryi istehlakı artımının demək olar ki, 1/5 hava şəraitinin dəyişməsi, yəni isti yay 
günlərində soyutma, soyuq hava şəratində iso əlavə istilik ehtiyacları ilə əlaqəli olmuşdur. 
Azərbaycan EES 21 iyul 2021-ci il tarixində ölkə üzrə havanın hərarətinin yüksək olması səbəbindən 
maksimum gücə olan tələbat və sutkalıq elektrik enerjisi istehsalında orta hesabla 12% artım qeydə 
alınmılşdır. 

Azərbaycan Dünyanın ən qədim neft və qaz hasil edən ölkələrindən biri olmaqla (2019-cu il 
üçün 1 milyard tondan çox təsdiqlənmiş neft ehtiyatına, həmçinin təxminən 2,1 trilyon kubmetr böyük 
təbil qaz ehtiyatlarına sahib olmaqla [3] yanaşı, bərpa olunan enerjinin böyük inkişaf potensialına 
malik olan ölkələr sırasına daxildir. Ölkə yüksək külək və günəş resurslarına və biokütlə, geotermal 
və su elektrik stansiyaları baxımdan əhəmiyyətli perspektivlərə malikdir -su elektrik stansiyaları üzrə 
5500 MVt, günəş FV >23000 MVt, külək (on-shor) x3000 MVt, tullantı / bioenerji 380 MVt, yekun 
olaraq 527000 MVt təşkil edir. Bundan əlavə, ucqar kənd ərazilərində müstəqil sistemlər üçün əlavə 
bərpa olunan istilik (biokütlə, günəş termal, geotermal) potensialı təqribən 0,6 ММЕТЛІ səviyyəsində 
qiymətləndirilir. Bütün ssenarilər üçün BEM-lərin əhəmiyyətli dərəcədə tətbiqi üçün proqnoz planlar 
hazırlanmışdır. Bu potensialın reallaşdırılması məqsədilə hökumət 2016-cı ildə qəbul etdiyi Strateji 
Yol Xəritəsində öz qarşısına 2020-cı ila qədər bərpa olunan enerji mənbəyinə 420 MVt potensialın 
əlavə edilməsini planlasdirimisdir. Bu məqsədi rəhbər tutaraq, BEM-in tikintisi üzrə bir neçə 
müqavilə imzalanmışdır (“Xızı” və “Pirəküşkül” KES, qoyuluş gücü 240 MVt, “Ələt” PV stansiyası, 
249 MVt) və Qarabağ bölgəsində qoyuluş gücü 300 MVt-da artıq KEQ və PVQ-m tikintisi üzrə ilkin 
texniki-əsaslandırma sənədləri icra olunmaqdadır. Eyni zamanda, elektrik enerjisinin 
hesablaşmasında “net metrinq" qaydasının tətbiqi yaxın 5 ildə fərdi evlərdə kiçik güclü DBEM əsaslı 
aktiv istehlakçıların meydana gəlməsini şərtləndirəcəkdir. Təqribi hesablamalara görə bu tip 
mənbələrin ümumi gücü 1000 MVt artıq ola bilər. Başqa sözlə, yaxın beş ildə DBEM-in йтйиті 
həcmi 2000 MVt-dan artıq proqnozlaşdırıla bilər ki, bu da EES-nın qoyluş gücündə payının 25%- 
dən çox təşkil edəcəyini proqnozlaşdırmağa imkan verir. 

Bərpa olunan enerji mənbələrinin enerji balansında daha yüksək payı iqtisadiyyatın 
şaxələndirilməsindən başqa, digər çoxsaylı faydalar gətirə və eləcə də BEM ənənəvi olaraq neft və 
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qazın üstünlük təşkil etdiyi ölkədə texnoloji innovasiyalara imkanlar yaratmaqla və iqtisadi dəyərin 
və bununla bağlı ÜDM artımının təmin edilməsi üçün yeni sektorlar açmaqla yeni məşğulluq 
imkanlarının katalizatoru qismində çıxış edə bilər. Bundan başqa, enerji effektivliyinin davamlı 
artması ilə yanaşı BEM-dən geniş istifadə ölkə daxilində neft və qaz istehlakının azalmasına gətirib 
çıxara bilər. 

Azərbaycanda enerjiyo tələbatın proqnozlaşdırılan artımını nəzərə alaraq, bərpa olunan enerji 
mənbələri ənənəvi mənbələrlə müqayisədə layihənin icrasının daha qısa müddətə başa gəlməsi 
səbəbindən enerji generasiya gücünü daha tez əldə etmək üçün texniki həll təmin edir. 

Bərpa olunan enerji mənbələri, həmçinin Azərbaycan Qlobal İqlim Dəyişmələri üzrə Çərçivə 
Konvensiyası üzrə götürdüyü hədəflərinə nail olmaq üçün kifayət qədər nəzərə çarpan aşağı karbonlu 
həll təklif edir. Ölkə Paris Sazişi çərçivəsində özünün Milli Səviyyədə Müəyyən Edilmiş Töhfə 
(NDC) sənədində qəbul edilmiş 1990-cı il baza ilinə nisbətən 2030-cu ilədək İQ emissiyasını 35% 
azaltmağı öhdəsinə götürmüşdür və bu hədəfə nail olmaq üçün BEM-dən istifadə xüsusilə vacib 
əhəmiyyət kəsb edir. 

Azərbaycan Respublikasının “Azəraycan 2030: sosial-qtisadi iksafinin Milli Prioritetlorindo 
“yaşıl artım” mərkəzi yerlərdən birini tutur. Eyn zamanda, Azərbaycanın da tərəfdaş olduğu Birləşmiş 
Millətlər Təşkilatının “Dünyamızın transformasiyası: 2030-cu ilədək dayanıqlı inkişaf sahəsində 
Gündəlik”dən irəli gələn öhdəliklərin icrası istiqamətində də xüsusi əhəmiyyət kəsb edir. 

Hazırda Azərbaycan EES-nin eneji mənbələrinin qoylus gücündə payına görə ənənəvi mənbələr 
tərbən 84%, BEM iso 16% təşil edir (şəki/ 7).O cümlədən, ənənəvi enerji mənbələri əsasən İES (40%), 
QTES (46%) və daxılıyanma mühərriklərindən (14%) ibarətdir. ВЕМ-ш struturunda SES əsas paya 


Enerji mənbələrinin strukturu, MVt /% BEM-nin mənbələrinin strukturu, MVt 


1318,93/ 
16,376 


6282,50/ 1180,38 


83,7% 


ш Ənənəvi ESt E Su elektrik st. = Külək elektrik st. m günəş elektrik st. Е Віоепегј 





Şək. 1. 01.01.2020-ci il tarixə (a) enerji mənələrini strukturu, (b) BEM-in struktru 


(90%) malikdir. Bərpa olunan enerji mənşəli generasiya güclərinin orta hesabla 90%-ini 
tonzimlononon BEM sinifinə daxil olan su elektrik stansiyaları təşkil edir (şəkil 1). 

Qeyd olunmalıdır ki, ölkə üzrə generasiya mənbələrinin qoyuluş gücündə “Azərenerji” ASC- 
nin payı 6,5 QV vo ya 86% olmaqla dövlət elekroenergetika sisteminin əsas enerji istehalçısı vo 
sistem operatoru funksiyalarını yerinə yetirir. Bu baxımdan, Dövlət elektroenergetika sisteminin 
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dayanaqlı idarə olunması ASC-yə həvalə olunmuşdur vo ВЕМ-ш hansı həcmdə və EES-in hansı 
nöqtələrinə və hansı formada inteqrə olunması üçün siyasəti məhz “Атогепегг ASC 
formalaşdırmalıdır. 


2. BEM-in EES-nə artan penetrasiya səviyyəsinin sistemə təsiri 


Elektroenergetika sistemlərinin normal parametrlər çərçivəsində dayanıqlı fəaliyyətinin təmin 
olunmasını xarakterizə edən vacib indikatorlardan biri enerjisistemlorinin çevikliyi hesab olunur. 
EES-də çeviklik, təklif-tələb balansını təmin etmək, gözlənilməz vəziyyətlərdə dayanıqlığı qorumaq 
və tələb-təklif tərəfindəki qeyri-müəyyənliyi aradan qaldırmaq qabiliyyətidir. Cevikliyi təmin etmək 
üçün yeni metod və idarəetmə tələbləri, enerji qeyri-müoyyonliyi və mövcudluğu ilo xarakterizə 
olunan bərpa olunan enerji ehtiyatlarının enerji balansında artması səbəbindən ortaya çıxmışdır. Bu 
təhlildə, EES- çevikliyi, çevik EES-nin xüsusiyyətləri, çevik mənbələri və qiymətləndirmə 
parametrləri beynəlxalq ədəbiyyatın bir hissəsi olaraq təqdim edilir. Həmçinin, dəyişən bərpa olunan 
enerji mənbələrinin nüfuzunun (penetrasiya səviyyəsinin) enerji sisteminin keçid prosesinə, statik, 
dinamik və tezlik dayanıqlığına təsiri keyfiyyətcə qiymətləndirilir, növbəti mərhələdə Azərbaycan 
EES-nin nümunəsində kəmiyyətcə qiymətləndirmək üçün təməlin yaradılmasına cəhdlər edilir, eləcə 
də çevikliyin təmin edilməsi üsulları qiymətləndirilir. 

İlk növbədə DBEM-ın ənənəvi mənbələrdən fərqləndirən beş əsas xüsusiyyətini qeyd edək: 

" DBEG-nın maksimum gücü real-zaman kəsiyində mövcud olan küləyin surəti və günəş 

işığına uyğun olraq dəyişir, 

= güclərin dəyişməsinin yalnız bir-neçə gün, təkmil proqnoz sistemləri olan hallarda iso bir- 

neçə saat əvvələ qiymətləndirilməsi mümkündür, 

= sistemə konvertor texnologiyalarının köməyi ilə qoşulur - bu hər-hansı generasiya 

mənbəyinin gözlənilməz açılmasından EES-nin dayanıqlığının təmin olunmasında 
əhəmiyyətli ola bilər, 

" ПОВЕС modul tiplidir və daha çox paylanmış formada sistemə tətbiq edillir; 

= tobi yanacaqdan fərqli olaraq, küləyin surəti və günəş şüası, nəql oluna bilmir və effektiv 

mənbələr daha çox yük mərkəzlərindən uzaq məsafələrdə yerləşir. 

DBEM-ın ümumi oxşarlıqlarının olmasına baxmayaraq, külək və günəş FV fərqli cəhətləri 
kifayət qədərdir. 

Global miqyasda BEA tərəfindən tövsiyə xarakterli və uzun müddətli planlaşdırmada istifadə 
üçün nəzərdə tutlmuş DBEM-nın penetrasiya səviyyiyəsinə görə təsnifləşdirilmiş bes Fazanı qeyd 
etmək olar və Azərbaycan üçün də faydalı ola bilər [4]: 


1. DBEM-in penetrasiya səviyyəsinin birinci fazası (<3-5%) kifayət qədər sadədir: DBEM-in 
bu səviyyəsində sistemə heç bir nəzərə çarpan təsiri yoxdur. Külək və ya günəş 
qurğularının ilk qurğulardan daha çox böyük bir sistemə quraşdırıldığı yerlərdə, onların 
çıxışı və dəyişkənliyi diqqətdən kənarda qalır. 

2. İkinci Fazada (5 — 15%)DBEM-nin təsiri nəzərə çarpır, lakin bəzi əməliyyat təcrübələrini 
təkmilləşdirməklə DBEM-in bu səviyyədə həcmi olduqca asanlıqla sistemə inteqrasiya 
oluna bilər. Məsələn, çevik elektrik stansiyalarının və tələbdə baş verən dəyişkənliyi 
səmərəli şəkildə planlaşdırmağa imkanların yaradılması üçün ОВЕМ-п isehsal gücünü 
proqnozlaşdırmaq üçün təkmil proqnozlaşdırma sistemi qurmaq lazım ola bilər. 

3. Üçüncü Fazada (15 — 2090), çevikliyin təmin olunması piroritet məsələ olaraq gündəmə 
gəlir. Bir sözlə, DBEM dəyişkənliyinin təsiri həm ümumi sistemin işində, həm də digər 
elektrik stansiyalarında hiss olunur. DBEM-in penetrasiya səviyyəsinin bu mərhələsində 
EES-ın çevikliyi məsələsi aktuallaşır. Bu kontekstdə çeviklik sistemin tələb-təklif 
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balansında yüksək dərəcədə qeyri-müəyyənlik vo dəyişkənliklə məşğul olma qabiliyyəti 
kımı başa düşülür. 

4. Dördüncü Fazada (20 -60%) yeni problemlər yaranır. Bu Fazada sistemin dayanıqlığının 
təmin olunması mühüm problem kimi meydana çıxır - çox qısa zaman kəsiyində 
həyacanlandırıcı təsirlərə sistemin dözmə qabiliyyəti. 

Sistem dayanıqlığlnın DBEM-ın inteqrasiyası üçün bir problem halına gəldiyi nöqtə bir-çox 
sistemlərə xasdır və onilliklər əvvəl qəbul edilmiş mühəndislik qərarlarından asılı ola bilər. Bu 
səbəbdən dayanıqlıq problemlərinin hansı penetrasiya səviyəsində və nə vaxt ortaya çıxacağını 
söyləmək üçün sadə bir qayda yoxdur. 

Enerji sistemi yuxarıda qeyd olunan bir mərhələdən digərinə dərhal keçməyəcək. Fazalar 
konseptualdır, sadəcə vəzifələrin prioritetlosdirilmosino kömək etmək məqsədi daşıyır. Məsələn, 
çevikliklə əlaqəli problemlər Üçüncü Faza halına gəlməzdən əvvəl, Ikinci Fazada tədricən ortaya 
çıxacaq. | 

Eynilə, sistem dayanıqlığı ilo bağlı bəzi problemlər artıq Üçüncü Fazada ortaya çıxacaq. 

Dördüncü Fazadan sonrakı beşinci və altıncı fazalar mövcud olsa da, hazırkı inkişaf 
mərhələsində baxılmır. 

Beləliklə, yuxarıda qeyd olunaları nəzərə alaraq, dəyişən bərpa olunan əsaslı generasiya 
mənbələrinin EES-ə penetrasiya səviyyələrinə görə təsiri necədir ? Hansı penetrasiya səviyyəsi 
optimal hesab oluna bilər? Bu və digər suallara sadə cavab yoxdur - iki eyni dinamik xassəli EES 
mövcud olmadığı kimi, iki fərqli ölkənin nə günəş, nə də külək mənbələri eyni xüsusiyyətlərə malik 
deyil, Nəticədə, külək və günəş enerjisinin müəyyən bir penetrasiya saviyasini və ya dəyəri ilo əlaqəli 
inteqral parameter və sadə qaydalar əldə etmək mümkün deyil. Eyni ilə, iki müxtəlif ölkə DBEM-in 
müxtəlif Fazlarında eyni BEM əsaslı generasiya payına malik ola bilər. Şərti olaraq DBEM-in 
penetrasiya səviyyələrinə görə Fazalara bölünməsi planlaşdırmada mümkün inteqrasiya səviyyəsini 
qiymətləndirmək üçün istifadə oluna bilər. Azərbaycan EES üzrə penetrasiya səviyyəsinin kritik 
həddinin ( təqribən 22096) təyin olunması, onun xüsusiyyətlərini nəzərə alamaqla, hər bir konkret 
layihə üzrə müvafiq hesablamaların aparılmasından sonra müəyyən oluna bilər. 

Eyni zamanda, aydındır ki, yüklərin dəyişən təbiətinə əlavə olunaraq BEM əsaslı generasiyanin 
anı, yüksək dəyişkənlik xüsusiyyəti ənənəvi güc mənbələrinin normal fəliyyətini kifayət qədər 
çətinləşdirir. 

Dəyişkənlik həm palnlaşdırma, hom də idarəolunma mərhələsində nəzərə alınmalıdır. 
Baxmayaraq ki, ənənəvi tələbatçıların istehlak gücü daima dəyişir, lakin uzun müddətli resurs 
planlaşmasında belə dəyişmə mühüm əhəmiyyətə malik deyil, DBEM-nin gündəlik tsikldə 
dəyişkənliyi iso bir gün erkən planlaşdırmada vacib faktor kimi çıxış edir. Cox qısa bir müddət olaraq 
təyin edilə bilən zaman aralığında (ms), DBEM əsaslı generasiyada anı dəyişikliklər səbəbindən yeni 
idarəetmə sistemləri tələb olunur. Bunlar aktiv və reaktiv güc, gərginlik və tezliyin dəyişmə dərəcəsi 
üçün idarəetmə sistemləridir [5]. 

DBEG-nın geniş tətbiqi halında əməliyyat çevikliyinin növü vaxt miqyasından asılıdır. Saniyə 
- dəqiqə vaxt aralığında daha çox tezlik nəzarəti və güc ehtiyatlarına ehtiyac olsa da, bir neçə 
dəqiqədən bir saata qədər və bir saatdan bir günə qədər vaxt aralığında erkən planlaşdırmanın 
çevikliyi tələb olunur. 

DBEM-də KQ və FVQ-ın üstünlük təşkil etdiyini nəzərə alaraq, ayrılıqda bu mənbələrin EES- 
ə təsirlərini nəzərədən keçirək. 
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2.1. Külək elektrik stansiyasının (KES) EES-in dayanıqlığına təsiri 


Holttinen'in araşdırmalarına görə, Skandinaviya ölkələri üçün külək enerjisinin EES-a 10% 
nüfuz etməsi, ehtiyat gücün həcminin quraşdırılmış külək elektrik stansiyalarının gücünün 1,570 -4% 
səviyyəsində artırılmasını tələb edir.Diggər tərəfdən, elektrik enerjisino olan tələbatı təmin edən 
elektrik stansiyalarının minimum buraxıla bilən gücü nə qədər aşağı olarsa, elektrik stansiyalarını 
dayandırmağa ehtiyac olmadan sistemə daha çox DBEG-nin nüfuz etməsi təmin edilə bilər. Buna 
görə də, artan güc ehtiyatı tələbləri ilə yanaşı, elektrik stansiyalarının buraxıla bilən minimum gücləri 
da qısa müddətli planlaşmada problemlərə səbəb ola bilər [6]. BEM-in penetrasiya səviyyəsi (5-35)% 
və müxtəlif dislokasiya nöqtələri üçün gərginliyin dayanıqlığı baxımından Naser və digərləri [7,8] 
tərəfindən tədqiqat işləri aparılmışdır. Müəyyən olunmuşdur ki, daha aşağı penetrasiya səviyəsində 
EES daha dayanıqlı fəaliyyət göstərir. Nəticələr göstərir ki, KES -in bir neçə nöqtədə şəbəkəyə 
qoşulması tək nöqtəli əlaqə ilə müqayisədə gərginliyin sabitliyinə müsbət təsir göstərir. Hossain və 
digərlərinə görə 191, İkiqat Qidalanan İnduksion Generator (Doubly Fed Induction Generator) (DFİG) 
turbinləri, Sinxron Generatorlara (SG) qədər reaktiv güc verə bilmir və böyük qısaqapanma 
cərəyanları yarada bılmır. Beləliklə, zəruri hallarda reaktiv güc ilə təmin edilən gərginlik dəstəyi 
DFIG üçün SG-don daha pisdir. Digər tərəfdən, Londero və digərləri [10], yüksək penetrasiya 
səviyyələrinin sistemin gərginlik üzrə dayanıqlığına müsbət töhfə verdiyi ifadə edilmişdir. 
Türbinlərin penetrasiya səviyyəsi artdıqca sistemə daha çox reaktiv güc dəstəyi verə biləcəyi ifadə 
edilir. 

Meegahapola vo Flynn [11], 39 düyünlü test sistemində çox yüksək (40%) külək enerjisinin 
(ОНО) nüfuz etməsinin enerji sisteminin keçid və tezlik dayanıqlığına təsirini araşdırmışdır. Küləyin 
yüksək nüfuz etdiyi bir bölgəyə yaxın ərazidə keçid prosesinin dayanıqlığına mənfi təsir etdiyi 
müşahidə edilmişdir. Bunun əsas səbəbi, KQ -lər qismən yükləndiyi zaman aktiv enerji istehsalının 
azalması və reaktiv güc istehlakının artması kimi göstərilmişdir. SG -lərə yaxın yerdə qəza halında, 
KQ-lorin sinxronlaşdırıcı gücünün artırılması yolu ilə şəbəkədəki güc axınlarına töhfə verməklə keçid 
prosesinin dayanıqlığının yaxşılaşdırdığı müşahidə edilmişdir. Edrah və başqaları 1121, bir SG və 
DFIG turbinlərindən ibarət külək stansiyasından istifadə edərək üç generatorlu 9 düyünlü test 
sisteminə müxtəlif ssenarilər tətbiq etməklə əldə edilən nəticələrə görə, sistemin rotor bucağının 
dayanıqlığına SG əvəzinə ekvivalent güc kimi DFIG turbinlorin istifadəsi mənfi təsir: göstərir. 
Bununla yanaşı, bu təsirin lazımı idarəetmə strategiyalarından istifadə etməklə azalda biləcəyi 
vurğulanır. 

Qureshi və İqbal [13], Rororu Qısaqapanmış İnduksion Generator (Squirrel Cage Induction 
Generator) (SCIG) vo DFIG turbinlorinin 9 düyün 3 masinli sistemindon istifado edorok sistemin 
tezlik üzrə dayanıqlığına təsiri üzrə araşdırmanın nəticələrinə görə, KQ-lorin tezlik meyletmosinin 
söndürülməsinə töhfəsi SG-lərə nisbətən SCİG-dən daha azdır. DFIG turbinlər, rotorun mexaniki 
sürəti şəbəkə tezliyindən ayrıldığı üçün tezlik meyletmosino cavab vermir. 


2.1.1. Külək elektrik stansiyaları interqrə olunmuş EES-do gərginliyə görə dayanıqlıq 
məsələləri 


Enerjisistemdo külək enerji mənbələrinin qoyuluş güclərinin artması ilo reaktiv güc və 
gərginlik üzrə rejimin tənzimlənmə effektivliyinin yüksəldilməsi yolu ilə gərginliyin dayanıqlıq 
səviyyəsinin artırılmasına nail oluna bilər. Sistemə qoşulmuş külək stansiyalarının güc hasilatının 
qeyri-müoyyonliyi və davamlı dəyişkənliyini nəzərə aldıqda sistemin dayanıqlı vəziyyətini təyin edən 
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parametrlər — şəbəkə üzrə ötürülüən güc həddi, onun düyünlərində gərginliyin kritik qiymətləri 
təsadüfi dəyişmə xarakterinə malik olacaq. Bu parametrlərinin qiymətlərinin təyin edilməsi üçün 
sistemin rejimlorinin ehtimal modelləşdirilməsi yolu ilə ehtimal güc paylanmasının (EGP) 
hesabatlarını aparmaq tələb olunur. Müasir dövrdə enerjisistemdo stoxastik yük və ya generasiya ilo 
şərtlənən EGP-nın müxtəlif ədədi və analitik analizi üsulları təklif edilmişdir [14]. 

Külək stansiyasının gərginliyin dayanıqlığına təsirinin tədqiqi  enerjisistemo inteqrasiya 
olunan KEQ-in qoyuluş güclərinin hədd qiymətlərinin təyini üçün, həm də bu tədqiqatların nəticələri 
əsasında istismar şəraitlərində sistemin fəaliyyətinin etibatlılığının saxlanması, başqa sözlə, KEQ-dən 
olan hasilatın fluktasiyaları zamanı elektrik sisteminin bütün düyünlərində gərginliyin dayanıqlığının 
saxlanması üçün vacibdir. КЕО-ш generasiya gücünün stoxastik dəyişmələri şəraitində sistemin 
rejiminin ehtimal modeli aşağıdakı kimi ifadə oluna bilər [14]: 
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Proq „0, KEQ( — КЕО-ш t zamanı ərzində aktiv və reaktiv güc hasılatları, 


P У ( E Bic - 1 düyünündə aktiv və reaktiv yüklər, 
Q SKB(i) - 1 düyününo qosulmus kompensasiya qurgularinin gücü; 


С,» В, = lvə J düyünlori arasindaki aktiv vo reaktic keciriciliklor; 


A — şəbəkənin yüklənmə əmsalı (yük artımı omsali); 


N - şəbəkə düyünlərinin ümumi sayı N = Npy + Npg ; 
U i U j- lvo J düyünlərindəki hesabı gərginliklər, 


щот zamanın k-c1 fiksə olunmuş qiyməti. 

Yuxarıda qeyd olunan ehtimal modelin köməyi ilo gərginliyin dayanıqlığı çoxsaylı təkrar 
hesabatlar ilə imitasiya olunan EGP nəticələrinin analizi əsasında təyin edilir. Aparılan hesabi 
eksperimentlərin alqoritmi 1141-də verilmişdir. 


2.1.2. Gərginliyin dayanıqlığı üzrə hesablam-təcrübə nəticələri 


Real Azərbaycan enerjisisteminin Yaşma düyünü üçün hesabi-tocrübi nəticələr əsasında güc 
hasilatının ehtimal xarakteristikaları, onun gərginliyin dayanıqlığına təsiri və həmçinin gərginliyin 
dayanıqlığının saxlanması şərtini təmin edən yekun külək parkının gücü təyin edilmişdir. 

Hesabat modelləşməsi külək sürətinin dəqiqəlik verilənləri əsasında aparılmışdır. Şəkil 2 və 3- 
də küləyin sürətinin və KEM-in hasilat gücünün dəyişmə fraqmentləri verilmişdir. Ölçmə periodu 
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sohor yük maksimumu dóvrüno uygundur. Diaqramdan góründüyü kimi verilmis periodda 
stansiyanin hasilatinin orta qiyməti 20 MVt təşkil edir (şəkil 3). 
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Sakil 3. KEQ-in giiciiniin dayisma diaqrami 


Enerjisistemin düyününo gosulan külok stansiyasinin maksimal mümkün gücünün qiymoti 
külək turbinlərinin müxtəlif yekun qiymətləri vo küləyin sürətinin stoxastik dəyişmələrinə uyğun 
EGP-nın analizi əsasında axtarılır. Şəkil 4-də nəzarət olunan elektrik veriliş xəttinin maksimal 
yüklənməsinin ümumiləşdirilmiş hesabat nəticələri və güc həddinin bu qiymətinə uyğun kritik 
gərginliyin qiyməti verilir. 

Şəkil 4-dən göründüyü kimi, baxılan sxem fraqmentinds gərginliyin dayanıqlığı təmin 
olunmaqla enerjisistemo qoşula biləcək KEM-in maksimal gücü 200 MVt və uyğun olaraq kritik 
gərginliyin qiyməti 64,496 təşkil edir. Müəyyən dayanıqlıq ehtiyatını nəzərə aldıqda bu gücü 50 MVt 
aşağı götürməklə, 150 MVt təyin etmək olar. Bu halda müəyyən dayanıqlıq ehtiyatı ilə baxılan halda 
gərginliyin dayanıqlığı təmin ediləcəkdir. Qeyd olunan metodika və alqoritm əsasında enerjisistemin 
digər prioritet düyünləri üçün analofi tədqiqatlar aparmaq olar. 
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2.2. Fotovoltaik sistemlərin EES-nə təsiri 


Kaliforniya Müstəqil Sistem Operatoru (CAISO) 1141 tərəfindən müxtəlif bilrləşmə sxemləri 
üçün FV sistemlərin gərginliyi və reaktiv güc reaksiyaları ilə bağlı araşdırmaya görə ötürmə 
şəbəkəsinə qoşulan FV-nın payının yüksək olması səbəbindən həddindən artıq ifrat gərginlik 
problemlərinin qaçılmaz olduğu ifadə edilmişdir. FV ilə bir sistemdə fəaliyyət göstərən STATCOM- 
ların ifrat gərginliyə səbəb olduğu da qeyd edilmişdir. Bunun əsas səbəbi, pozuntunun aradan 
qaldırılmasından sonra aşağı işləmə sürətinə malik olan STATCOM-lar tərəfindən sistemə bir neçə 
dövr ərzində reaktiv gücün verilməsidir. Eftekharneyad və digərlərinə 1151 görə ən çox təsirlənən 
sistem parametri gərginliyin amplitududur. Xüsusilo,2046 va daha yüksək penetrasiya səviyyələrində 
ötürmə xəttin birləşmə şinlərində həddindən artıq gərginlik meyletmosi mydana gəlir. Həmçinin 
keçid prosesi müddətində gərginlik uçqunu baş verir. 

Tamimi və başqaları, mərkəzləşdirilmiş və paylanmış PV sistemlərinin Ontario enerji 
sisteminin gərginlik üzrə dayanıqlığa təsirlərinin araşdırmasının nəticələrinə görə (quraşdırılmış 
gücü 2000 MVt -a qədər olan müxtəlif penetrasiya səviyyələri araşdırılmışdır), paylanmış FV -nin 
mərkəzləşdirilmiş sistemlərlə müqayisədə gərginliyin dayanıqlığını əhəmiyyətli dərəcədə artıra 
biləcəyini müəyyən edilmişdir [16]. 

Keçid prosesinin dayanıqlığının yoxlanılması məqsədilə aparılan todqiqatlrda, penetrasiya 
səviyyəsinin təsiri ilk olaraq 5-30% aralığında araşdırılmalardan əldə olunan nəticələrə görə, FV-nin 
10%-don yuxarı penetrasiya səviyyələri üçün sistemin dayanıqlığının yaxşılaşdırılması bu elektrik 
stansiyalarının qısaqapanma müddətində əsas rejim parametrlərinin tənzimlənməsində töhfə vermə 
imkanlarından asılıdır [14]. Eftekharnejad və digərləri 1151 fərdi mənzillərdə quraşdırılan kiçik və 
sənaye tipli FV sistemlərinin penetrasiya səviyyəsinin bir-biri ilə əlaqəli böyük bir enerji sisteminə 
təsirini araşdırmışdır. Bu məqsədlə, ənənəvi generasiya payı 50% -ə qədər azaldılaraq, 
FV penetrasiya səviyyələri araşdırılmışdır. Təhlillər göstərmişdir ki, yüksək FV penetrasiya səviyyəsi 
sistemin keçid prosesinin dayanıqlığına müsbət və mənfi təsir göstərir. Həmçinin,FV penetrasiya 
səviyyələri, sistem topologiyası, pozuntunun növü və yeri təsirin xarakteri üçün əsas amillərdir 
(müsbət və ya mənfi). Bundan əlavə, fərdi evlərdə FV -nin böyük bir hissəsinin fəaliyyətinin 
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dayanması yaxınlıqdakı SG-larin rotor bucaqlarının meyletməsinə və gərginlik fluktasiyalarına səbəb 
olmuşdur. 

Eftekharneyad və digərləri öz tədqiqatlarında 69 КУ gərginlik səviyyəsində toplanmış böyük 
ölçülü FV və fərdi evlərdə quraşdırılmış FV-lərdən istifadə etmişdir. Əldə olunan nəticələrə görə, 
böyük şərti ənənəvi generatorların FV mənbələri ilo penetrasiya səviyyəsi 30% -dən 40% -ə qədər 
artdıqca dempfer nisbətində əhəmiyyətli azalma müşahidə olunmuşdur. Penetrasiya səviyyəsinin 
artması, sistemin inersiyasinin azalmasına səbəb olaraq sistemin dempferləmə qabiliyyətinin 
azalmasına səbəb olmuşdur 1171. 
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Şək.2. Azərbayca EES-nin inersiya sabitinin dinamikası 


Şəkil 2-dən aydın görsənir ki, 2020-2030-cu illərdə Azərbaycan ЕЕЅ-аә BEM-in geniş 
penetrasiya səviyyəsi ilo əlaqəli ekvivalent inersiya göstəricisi hesabat dövrünün sonunda 18,4% 
azamağa məruz qalır. Nəticədə, EES-də güc balansının gözlənilməz pozulması halında inersiya 
reaksiyası azalır və belə vəziyyət sistemin çevikliyi və enerji təhlükəsizliyi ilə bağlı yeni problemləri 
gündəliyə gətirir. Gündəlik kontekstdə nəzərdən keçirsək, 21 -ci əsrin sisteminin ətaləti, gecə 
satlarında aşağıdır (tələbat az olduqda, az elektrik stansiyası işləyir) və yüksək tələbat zamanı 
yüksəkdir. FV-nin sistemin ətalətinə təsiri, həm yüksək, həm də aşağı tələbat dövründə baş verəcək 
küləkdən daha az olacaqdır. 

Aydındır ki, enerji balansında BEM-in payının artması, sistemin çevikliyinə təsirini 
qiymətləndirmək üçün təcili bir ehtiyac yaradır, çünki günəş və külək üzərində qurulan generasiya 
"kəsilən" və gözlənilməz bir təbiətə malikdir. Tezliyin hər iki (birinci və ikinci) tənzimlənmə sistemi 
inersiya reaksiyası По əlaqəli olduğu üçün ВЕМ-ш inteqrasiya səviyyəsindən ciddi şəkildə asılıdır. 

Alquthamı və digərləri [18] FV penetrasiya səviyyəsinin səviyyələrini 5%, 1096vo 20% hədləri 
üçün tezlik üzrə dayanıqlıqa bağlı digər araşdırmalarında göstərmişdir ki, FV-nin 2046 penetrasiya 
səviyyəsindən başlayaraq EES-nin tezliyə görə dayanıqlıq qabiliyyəti kifayət qədər azalmağa məruz 
qalır. 


Nəticə 


1. EHS-də DBEG-nın penetrasiya səviyyəsinin artması (72096) sistemin davranışına birbaşa 
təsir göstərir. DBEG-nin buraxılış gücündəki rəqsləri EES-nin statik, dinamik 


dayanıqlığına, keçid prosesinə, tezlik və gərginliyin tənzimlənməsinə mənfi təsir göstərə 
bilər. 


10. 


11. 


12. 
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Fazalari ümumilosdirorok, bóyük bir enerji sisteminin müoyyon hissolorinin sistemin qalan 
hissəsindən əvvəl daha inkişaf etmiş bir Fazaya girə biləcəyi mükündür. Ənənəvi olaraq, alt 
şəbəkələrin bir-biri ilə əlaqəli əsas şəbəkəyə qoşulduğu bir vəziyyətdədir. 

DBEM çıxışının vaxtının enerji tələbi ilo korelyasiyasi,sistemin ölçüsü və sistem daxili 
əlaqələri və operatorun DBEM çıxışını proqnozlaşdırma qabiliyyəti, bütün bunlar yeni bir 
inteqrasiya fazasının nə vaxt baş verəcəyini müəyyən edəcək. 

DBEM-nın səmərəli şəbəkə inteqrasiyası, yerləşmə Fazasına uyğun və mütənasib olan 
tədbirləri üzə çıxaracaqdır. Bəzi sistemlərdə bu tədbirlər mövcud aktivlərdən istifadə 
etməklə həyata keçirilə bilər, digərlərində isə əlavə infrastruktura investisiya qoyulması 
tələb oluna bilər. Hər iki halda, artan DBEM payına uyğunlaşmaq uzun müddətdə daha 
böyük xərclərə səbəb olacaq və enerji sisteminin təhlükəsizliyini təhdid edə bilər. Əksinə, 
həddindən artıq yüksək tələblərin qoyulması xərcləri və/və ya DBEM-ın tətbiq sürətini 
azalda bilər. 

Enerji təhlükəsizliyi və çevikliyi üçün kritik olan və sistemin səmərəli fəaliyyəti üçün lazım 
olan tədbirləri ayırmaq məqsədəuyğun hesab olunur. 

İkinci Fazadan başlayaraq sistem operatorunun DBEM generasiyanın böyük bir hissəsini 
dayandıra bilməsi üçün səlahiyyətin verilməsi vacibdir. Əks təqdirdə, enerji təhlükəsizliyi 
risk altında ola bilər. Əksinə, təsirli bir proqnozlaşdırma sisteminin olmaması sistem 
operatorunu digər generatorları həddindən artıq yükləməyə və sonra külək istehsalını 
həddindən artıq azaltmağa-xərcləri artırmağa səbəb ola bilər. 

Fərdi evlərdə quraşdırılmış günəş invertorlarinin tənzimləmə parametrlərin dəyişdirilməsi 
imkannin olması zəruridir. 

Enerji saxlanclarının tətbiqinin diqqətdə saxlanılması vacib hesab olunur. 

Şəbəkələri balanslaşdırmaq üçün DBEM-dən istifadə olunması imkanlarının araşdırlması 
tələb olunur. 

Sistemin tənzimləmə qabiliyyətinin artırılması məqsədilə, “Tovuz” və digər böyük SES-in, 
eləcə də Hidroakumlasiya stansiyasının tikintisinin planlaşdırılması məqbul hesab olunur . 
Ötən dövrdə, DBEM-lər enerji sektoruna sistem və köməkçi xidmətlər göstərmək 
öhdəliyindən “azad edilmişdi”. Artıq, sistemin etibarlılığını təmin etmək üçün xidmətlərin 
göstərilməsi və tolob-toklifin balanslaşdırılmasında iştirakının təmin olunması (məsələn, 
Danimarkada, bütün külək enerjisi qurğuları proqnozlaşdırılan gücünü təmin etməli vo 
çıxışı proqnozla uyğun gəlmirsə cərimələr ödəməlidir, bu da külək elektrik stansiyalarının 
sahiblərini proqnozlaşdırma vasitələrinin təkmilləşdirilməsinə sərmayə qoymağa vadar 
edir. Cənubi Avstraliyada tənzimləyici gələcəkdə tezlik nəzarəti üçün DBEM-dən istifadə 
edilməsinə qərar vermişdir). 

Ötürmə şəbəkəsinin gücləndirilməsi və genişləndirilməsi, paylanmış generasiyanın 
balanslaşdırılması, qonşu ölkələrlə qarşılıqlı əlaqələrin inkişaf etdirilməsi və beynəlxalq 
enerji bazarlarına çıxışın əldə olunması kimi tədbirlər, BEM-nin artan penetrasiya 
səviyyəsinin şəbəkəyə inteqrasiyası üçün çevikliyin artırılması ən yaxşı texniki həlldir. 
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DIAGNOSIS RESEARCH OF THE FAULT OF OIL-IMMERSED POWER TRANSFORMERS 
BASED ON PROBABILISTIC NEURAL NETWORK 


Atrafi Giyasi M.A., Bodurov V.C. Kerimov K.A. 


Abstract 


The oil-immersed power transformer 1s one of the key devices in the power system, they are filled 
with a fluid that serves various purposes. The fluid acts as a dielectric media, an insulator and a heat 
transfer agent. The most common type of fluid used in transformers is of a mineral oil origin. During 
normal operations, there is usually a slow degradation of the mineral oil to yield certain gases that are 
dissolved in the oil. However, when there 1s an electrical fault within a transformer, gases are generated 
at a much more rapid rate. DGA 1s probably the most widely employed preventative maintenance 
technique in use today to monitor on-line transformer operations, and a number of DGA interpretation 
guidelines have been developed by different organizations, e.g. IEC60559 [1], IEEE C57.104-1991 [2], 
CIGRE TF 15.01.01 [3] and GB7252- 87 [4]. The reliability of the power grid 1s guaranteed by its safe 
operation. Therefore, it is necessary to reduce transformer failures with precautionary measures. 
Nowadays, the three-ratio method of dissolved gas analysis (DGA) is the most effective and convenient 
method in transformer fault diagnosis. However, when using the three-ratio method as the interpretation, 
it exists some disadvantages such as coding defects and threshold criterion defect. A new way for this 
problem is provided by artificial neural network, which has the advantages such as parallel processing, 
self-adaptation self-study, association memory, non-linear mapping, and other features. The renewed 
thrust in smart power system networks along with the development of advanced methods in monitoring 
and diagnostics has resulted in major impetus to research in the related domain. 


The transformer internal faults are predicted in this paper by using a PNN algorithm. The high 
accuracy of diagnosis 1s confirmed by simulation results in the KNIME platform. 


Keywords: DGA, Two-ratio method, Rogers method, Diagnosis of gases dissolved in 
transformer oil, Power Transformer, Internal Faults in windings, Fault Diagnosis, Neural Network. 


ДИАГНОЗНОЕ ИССЛЕДОВАНИЕ НЕИСПРАВНОСТИ МАСЛОЗАГРУЖЕННЫХ 
СИЛОВЫХ ТРАНСФОРМАТОРОВ НА ОСНОВЕ ВЕРОЯТНОСТНОЙ НЕЙРОННОЙ 
СЕТИСЕТИ 


Атрафи Гияси М.А., Бодуров В.С. Керимов К.А. 
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Абстракт 


Масляный трансформатор - одно из ключевых устройств в энергосистеме, они заполнены 
жидкостью, которая служит различным целям. Жидкость действует как диэлектрическая среда, 
изолятор и теплоноситель. Чаще всего в трансформаторах используется жидкость на основе 
минерального масла. Во время нормальной работы обычно происходит медленное разложение 
минерального масла с выделением определенных газов, растворенных в масле. Однако, когда в 
трансформаторе возникает электрическая неисправность, газы выделяются гораздо быстрее. 
Анализ растворенного газа (DGA), вероятно, является наиболее широко используемым методом 
профилактического обслуживания, используемым сегодня для мониторинга работы 
трансформатора в режиме онлайн, и ряд руководств по интерпретации ОСА был разработан 
различными организациями, например IEC60559 [1], IEEE C57.104-1991 [2], CIGRE ТЕ 15.01.01 
[3] и GB7252-87 [4]. Надежность электросети гарантируется ee безопасной эксплуатацией. 
Следовательно, необходимо уменьшить количество отказов трансформатора с помощью мер 
предосторожности. В настоящее время трехуровневый метод анализа растворенных газов (DGA) 
является наиболее эффективным и удобным методом диагностики неисправностей 
трансформаторов. Однако при использовании метода трех соотношений в качестве 
интерпретации существуют некоторые недостатки, такие как дефекты кодирования и дефект 
порогового критерия. Новый способ решения этой проблемы предоставляется искусственной 
нейронной сетью, которая имеет такие преимущества, как параллельная обработка, 
самоадаптация, самообучение, ассоциативная память, нелинейное отображение и другие 
особенности. Возобновление популярности сетей интеллектуальных энергосистем наряду с 
развитием передовых методов мониторинга и диагностики привело к значительному 
стимулированию исследований в соответствующей области. 


В этой статье для прогнозирования внутренних неисправностей трансформатора 
используется вероятностный алгоритм нейронной сети. Высокая точность диагностики 
подтверждена результатами моделирования на платформе KNIME. 


Ключевые слова: ОСА, Двухступенчатый метод, Метод Роджерса, Диагностика газов, 
Растворенных в трансформаторном масле, Силовой трансформатор, Внутренние неисправности 
в обмотках, Диагностика неисправностей, Нейронная сеть. 


1. Introduction 


The transformer, which is the most important and fundamental component in a complex power 
system shall be responsible for power transmission and voltage converting. While the multi-protection 
outside and inside shield transformers, transformer failures still occur. Meanwhile, the frequency of 
accidents 1s quite high due to the design, manufacturing, technology, raw materials of defects, and natural 
disasters, etc. As a result, it is important to forecast and diagnose transformer faults in real-time. 
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Dissolved gas analysis (DGA) is the most effective and convenient method in transformer fault 
diagnosis. And there are three improved methods based on DGA technology, including the Two-ratio 
method proposed by Dornenburg (to distinguish heat and electricity) and the Rogers method (four-ratio 
method) [5,6]. But many questions in engineering practice are exposed, some of which are incomplete 
coding, absolute encoding boundaries, and the failure of faults diagnosis. As a result of the existence of 
these problems, a variety of artificial intelligence technologies for power transformer fault diagnosis have 
been proposed by domestic and foreign scholars. 


For example, the sample 1s to be classified effectively with the combination of genetic algorithm 
and FCM ) [7]. The optimization parameters of the genetic algorithm can be accurately found in a wide 
range by using SVM ) [8]. 


In the meantime, artificial intelligence algorithms such as expert system ) [9], fuzzy theory [10], 
rough set theory [11], evidence theory [12] are applied to the diagnosis of the power transformer. 
Nevertheless, because of the new fault produced by the power transformer in the actual process, the 
above methods can’t be updated in the training sample in real-time. PNN (PNN), a diagnosis of a network 
system with the abilities of strong fault tolerance and adaptive structure, to diagnose power transformers, 
has been adopted in this paper. 


The Three-ratio method is an effective interpretation method, but this method fails to clearly identify 
faults in transformers in borderline cases. At the same time, it 1$ inescapable to make a wide evaluation 
based on the test results of electrical test items, so as to make a correct interpretation of the nature and 
location of the fault. 


According to the different ratios of the volume fractions of gases, the given codes in Tab.1 


According to *DL/T722-2014 Guide to the analysis and the diagnosis of gases dissolved in transformer 
oll”, we can find out the transformer fault type and position. 


Table 1. 


Coding rules of three ratio method 
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2.Statement of Problem 


2.1 | Oil-immersed transformer fault diagnosis based оп РММ 
a. Pose problem 


The workflow was set up on the KNIME platform in order to effectively predict power 
transformer faults and the types of fault. The input sample data set D, which is expressed as (vi, v5,, Vn; 
К) is a previous fault data set. у; is the content of the characteristic gases in power transformers and Fis 
the corresponding fault types Fig.1. The input test data set X, which the types of fault and the accuracy 
of the prediction to be observed, only including the content of the gases and expressed as (vi, у2,, Vn). 


b. Data pre-processing 


The normalization processing, which is the linear function transformation method, is used 
in this article to make the data fall into the range of [0, 1] and to avoid the unequally using of the 
input variable with different physical meaning and dimension. It can also elimination the influence 
of dimensional, make the data become scalar, improve the accuracy of fault diagnosis and 
accelerate the convergence of training network. 


с. РММ algorithm process 


PNN is composed of three layers: the input layer, the hidden layer, and the output layer. This 
neural network model has the following features: the topology, connection weights and thresholds 
can be set immediately when training samples are attainable. Its training is concise and 
classification ability 1s strong. 


l. PNN diagnosis process 
e Connection weights matrix of the network УИ, W”° and training sample set D are read; 


ә Training sample set D is entered into the hidden layer and the output of the hidden layer 
is calculated; A nonlinear operator of hidden layer by gauss function as fellow: 


fk = ехр(-пеђ/о?) (1) 


o is smoothing factor. 


nety2 ŞI x, — wi (2) 
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smoothing factor o1 = oz =............. ON = 00. 


Probability sum O; ,which belongs to the first category,with the value as output from hidden 
layer to output layer is calculated; 


Ој= Xe (мкр. Hy) (3) 


J. Computation 


Failure prediction reasoning path of predicted data along the level of the network and the 
final output matching inference results. The output of failure data expressed as:(vi, və,, у»; 2), Vi is 
property values ,Z 1s predicting fault type. 





Fig. 1 PNN topological structure 


2.2 Simulation and analysis of Power transformer fault diagnosis based on PNN 


Dissolved gas analysis (DGA) is one of the most effective and important methods in oil- 
immersed power trans- former fault diagnosis. Among various diagnostic methods of DGA, the 
three-ratio method improved has the highest accuracy. The three-ratio of oil dissolved gas content 
are selected as input vectors of neural network. And the output vectors are the types of transformer 
faults. 
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2.3 Тһе sample data set 


A total of 300 sample data is got from the UCI machine learning repository, which contains seven 
attributes Tab.2 : 


CH4, C2H6, С На, С Но, Hə, CO, CO». 


Table 2 
Failure type and number: 
Failure type Number 
Low temperature overheat 52 
hyperthermia and superheating 260 
Low energy discharge 119 
high energy discharging 23 


Discharge and overheating 30 





From this 300 sample data, 250 sample data are selected randomly for PNN training, others are 
the test data for testing the accuracy of power transformer fault diagnosis which is based on PNN. 


2.4 Simulation workflow setting 


KNIME, a data implementation software, is used for simulation experiment and the 
simulation of fault diagnosis 1s realized by constructing workflows. 


The implementation steps as fellow: 


e 500 oil-immersed transformer fault data for database are read; 
e With the normalization processing, this data are divided into sample data and test data randomly; 
e PNN learning is used to learn of the data, and the fault prediction model for failure prediction 
is to be establish; the predict result and the existing result are compared, and the accuracy of 
prediction 1$ calculated. 
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3. Simulation results analysis 


With different algorithms, the sample data are learning and used to predict fault. simulation results 
as shown in Tab 3. 


From the above calculation results, the fault diagnosis accuracy based on PNN is 80%. It is higher 
than accuracy rate that based on other intelligent algorithms. As a result, satisfactory results on the 
accuracy of oil-immersed power transformer internal fault diagnosis are achieved by using РММ Fig.2 . 


PNN Learner (DDA) 


Normalizer 






Column Filter Partitioning 
Database Reader Color Manager — 







"9 œ 


pc Й spilt data 

read sample data — Fault distinguish P 
PNN Predictor Scorer 

Column Filter Row Filter 

Database Writer 


Fig. 2 Fault diagnosis workflow 


Table 3 


Fault algorithm and accuracy: 


Algorithm Accuracy prediction 
= = 7 
BP-neural network 0.767 


fuzzy algorithm 0.680 
m 
Bayesian algorithm 0.530 
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Conclusion: 


From the literature reviews, several fault diagnosis methods are required to detect all the faults 
of the Power Transformer. 


One of the methods based on PNN to establish workflow about oil-immersed power transformer internal 
fault diagnosis 1s introduced in this paper. With the ability of learning knowledge from the examples and 
convenient parameter settings, neural network with good fault-tolerance 1s obtained. The network, which 
has simple processing and rapid convergence, don't need training. Its corresponding weighting is the 
distribution of sample. The transformer fault sample library, which would change with the increase of 
transformer faults, should to be established. And with the full using of the strong ability of PNNs sample 
complementing, the practical application of diagnostic accuracy 15 increased. 
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DETERMINATION OF DENSITY DEPENDENCE ON TEMPERATURE OF GEOTHERMAL 
WATERS OF JARLI KURDAMIR REGION OF AZERBAIJAN 


AZORBAYCANIN KÜRDOMIR RAYONUNUN CARLI GEOTERMAL SUYUNUN 
SIXLIGININ TEMPERATURDAN ASILILIGININ TOYINI 
Babayeva S. Sh. DETERMINATION OF DENSITY DEPENDENCE ON TEMP... 


ОПРЕДЕЛЕНИЕ ЗАВИСИМОСТИ ПЛОТНОСТИ OT ТЕМПЕРАТУРЫ l'EOTEP- 
МАЛЬНЫХ ВОД ДЖАРЛЫ КЮРДАМИРСКОГО РАЙОНА АЗЕРБАЙДЖАНА 


Babayeva S. Sh. 
Azerbaijan Oil and Industry State University 


Annotation 

In the article, based on the analysis of literary sources, the dynamics of the development of 
alternative energy in different years is studied and a diagram 15 built for comparison based on the results. 
Based on the IRENA database, the dynamics of the growth of the potential of geothermal energy in the 
world countries has been determined. Explored the potential of geothermal energy in Azerbaijan. The 
density of a water sample taken from the Jarli geothermal spring in the Kurdamir region was 
experimentally investigated at atmospheric pressure and various temperatures (T = 278.15-343.15K), 
and the results were expressed by a polynomial equation. 

Keywords: alternative energy, Jarli, geothermal water, thermophysical characteristics, density, 

temperature dependence. 


Babayeva S.S. 
Azorbaycan Dóvlot Neft vo Sonaye Universiteti 


Annotasiya 

Məqalədə ədəbiyyatların analizinə əsasən, alternativ enerjinin illər üzrə inkişaf dinamikasi 
tədqiq edilmiş və nəticələrə əsasən müqayisə üçün diaqram qurulmuşdur. İRENA-nın məlumat bazasına 
əsaslanaraq, dünya dövlətlərində geotermal enerjinin potensialının inkişaf dinamikasının artımı təyin 
edilmişdir. Azərbaycanın geotermal enerji potensialı araşdırılmışdır. Kürdəmir”. rayonunun Carli 
geotermal mənbəyindən götürülən su nümunəsinin sıxlıqları atmosfer təzyiqində və müxtəlif 
temperaturlarda (T=278,15-343,15K) təcrübi tədqiq edilmiş və alınmış nəticələr polinomial tənlik 
vasitəsilə ifadə edilmişdir. 

Açar sözlər: Alternativ enerji, Carli, geotermal su, istilik-fiziki xarakteristikalar, sıxlıq, 

temperaturdan asılılıq. 
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Аннотация 

В статье на основе анализа литературных источников изучена динамика развития 
альтернативной энергетики в разные годы и построена диаграмма для сравнения по результатам. 
На основе базы данных IRENA определена динамика роста потенциала геотермальной энергии B 
странах мира. Изучен потенциал геотермальной энергии Азербайджана. Плотность пробы воды, 
взятой из геотермального источника Джарли в районе Кюрдамир, была экспериментально 
исследована при атмосферном давлении и различных температурах (T = 278,15-343,15К), и 
полученные результаты были выражены полиномиальным уравнением. 

Ключевые “слова: альтернативная энергетика, Джарлы, геотермальная вода, 

теплофизические характерис-тики, плотность, температурная зависимость. 


INTRODUCTION 


The depletion of fossil fuels and their negative environmental impacts during combustion have 
prompted many countries to use renewable energy sources. This, in turn, has led significant progress in 
the development of renewable energy technologies. [1-5]. 

The use of alternative energy in the world increased from 1.3% to 30% since 2001 till 
2020. Fossil fuel use declined from 85.2% to 20.1% over the corresponding years. [6]. 

The dynamics of growth in the total use of renewable energy in the world countries 1s shown in 
a diagram in Picture | based on the database presented [7] 
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Picture |. Total use of renewable energy in the world 
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Picture 2 shows the increase in the installed capacity of geothermal power plants in the world in 
recent years. This Picture is built based on the database shown in [7]. 
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Picture 2. Annual growth of installed capacity of geothermal power plants in the world 


Geothermal energy is an important part of the development of renewable energy. In 1904, for 
the first time in Italy, steam from a geothermal source was used to generate electricity in a steam power 
plant. Currently, the capacity of the Larderello Geothermal Power Plant in Italy is about 400 
MVt. Geothermal energy has developed in a number of countries. In recent years, the use of geothermal 
resources has been characterized by an increase in volumes and expansion of areas. Geothermal 
technologies dominate the energy balance of a number of countries, and the share of geothermal energy 
in the global energy balance is constantly growing. The leading position in the world is occupied by the 
USA (more than 4046 of the existing capacities). The installed capacity of the geothermal power plant in 
the Philippines is about 900 MVt, in Mexico - 700, in Italy - 500, in Russia - 21 MVt. [8]. In addition to 
generating electricity, geothermal energy can also be used for heating. The use of geothermal energy for 
the development of heating systems and the development of new schemes is on the agenda. 


MAIN PART 


Azerbaijan is also a country rich in mineral and thermal waters. Recently, interest in geothermal 
sources has been growing in Azerbaijan, as well as in all areas of using alternative energy. The potential 
of geothermal energy in Azerbaijan is estimated at about 800 MW. The area of application of these waters 
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is extensive and covers agriculture, medical, industrial and municipal heating. That is why scientists are 
very interested in studying geothermal waters as an object of research. In Azerbaijan, significant 
scientific work in the field of geothermal energy and mineral water deposits has been carried out for a 
long time [9-12]. Chemical analysis of their composition and etc. significantly studied. After the 
recognition of After the recognition of Azerbaijan's independence, research on the thermal and physical 
properties of geothermal energy water sources began. [13-14]. 

The Azerbaijan State Oil and Industry University, together with the University of Rostock in 
Germany, began the study of the thermophysical properties of geothermal energy sources as a new 
scientific direction. The main goal of this study is to identify the potential of alternative energy sources 
in the regions of Azerbaijan, including the Kura basin, and its use in the energy sector. 

The potential of geothermal energy sources located in the Kura basin in Azerbaijan has very 
promising parameters for energy use. The well, located in the village of Jarli, Kurdamir region, is 
considered the main geothermal source. This well attracts attention with its abundant water. The 
temperature of the discharged water is 95-98 °C, and the flow rate is 20,000 m? / day. [9] Water of Jarli 
is a layer of porous water. The chemical composition of this water has been studied by different 
researchers at different times. [6,15,16]. However, the thermophysical properties of this water have not 
yet been studied. 

We have investigated some thermophysical parameters of Jarli'S water experimentally. The 
density of geothermal water was measured using a DMA 5000 device. This device 1s one of the latest 
devices of the Anton Paar generation. The Anton Paar generation devices are fast, accurate, reliable, 
repeatable and have a modern user interface. 

The DMA 5000 high-precision instrument was used to determine the density and speed of sound 
of thermal water in the Jarli village of the Kurdamir region at atmospheric pressure. The experiments 
were carried out using a DMA 5000 densimeter (Picture 3) from Anton-Paar (Austria). The DMA 5000 
vibrating tube densimeter is based on the one-sided dependence of the oscillation period of the U-tube 
on the tube mass[17]. 





Picture 3. DMA 5000 device 
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One of the latest Anton Paar models, the DMA 5000 can simultaneously set 2 independent 
physical parameters, density and sound velocity with high precision. This device, which can measure 
two functions, has chambers for measuring density and speed of sound. Determines high-precision sound 
copy using the U-shaped tube method. The temperature in both chambers is controlled by a Peltie 
thermostat. 


1020 
1015 
1010 
1005 
1000 

993 


990 € 





Picture 4. Jarli geothermal water's dependence p - f (t) 


Picture 4 shows the temperature dependence of the density of Jarli's water. As the water 
temperature increased from 5 °C to 80 °C, its density decreased. 


p = -0,0037t^- 0,1165t+1019 (1) 


Table 1. Comparison of the density of the Jarli thermal water with the values obtained from the 
experiment and the values obtained from equation (1) 


°С Density of Jarli 
Density of Jarli — mE | 
thermal water, p, ош. P | Ao/pa -100% 
ƏTADƏTİTY1 
(təcrübə ) kq/m3 У 


kq/m’ 


5 1 1028,076 1018,525 9.551 
1017.905 1017.665 
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RESULT 





Studies have shown that the temperature dependence of the geothermal water density Jarli in 
Kurdamir region of Azerbaijan has the same dynamics as the temperature dependence of the density of 
fresh water and decreases rapidly with increasing temperature. The density of Jarli geothermal water is 
slightly higher than the density of fresh water. This is due to the fact that geothermal waters are rich in 
chemical minerals. According to the results of the experiments, the temperature dependence of the 
density of Jarli geothermal waters р / kg “ m-3 varies from Т = 298.15 К to high temperatures in the 
appropriate sequence, depending on the amount of chemical minerals in their composition. But in the T 
= (278.15-298.15) K range, this correspondence is violated. This discrepancy is associated with the low 
solubility of chemical minerals in geothermal waters at low temperatures compared to high temperatures, 
as well as the influence of the temperature dependence of the density of dissolved gases in the sample on 
water and mixtures with excess water content. 
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УДК 681.05:621 
Рустамов К.О. 


УМЕНЬШЕНИЕ ЭНЕРГОПОТРЕБЛЕНИЯ ОПТИМАЛЬНОЙ ПО 
БЫСТРОДЕЙСТВИЮ СИСТЕМЫ УПРАВЛЕНИЯ 


Азербайджанский Технический Университет 


В статье исследуется влияние вибрации управления на расход энергии в 
оптимальной по быстродействию системе и робастные свойства модифицированной 
системы управления. Решение классической задачи оптимального управления по 
критерию быстродействия приводит к двухпозиционному релейному управлению. В 
этом случае имеется следующие недостатки: в малой окрестности положения равновесия 
происходит вибрация управляющего сигнала, что приводит к перерасходу энергии; 
функция переключения зависит от времени и не имеет обратной связи. Учитывая эти 
недостатки предложены различные модификации постановок задач. Эти постановки 
отличаются многокритериальностью и добавлением функциональных ограничений. 
«Исправление» исходной задачи таким образом приводят к  усложнению 
математического аппарата и системотехнической реализации системы. Одним из 
конструктивных подходов, который использован в данной работе, является системный 
подход. При этом осуществляется линеаризация оптимального алгоритма: релейное 
управление заменяется на линейно-релейный, который в малой окрестности положения 
равновесия позволяет подключить обратную связь; нелинейная функция переключения 
аппроксимируется гиперплоскостью, в уравнению которого входят переменные 
состояния, кроме того упрощается задача синтеза для систем высокого порядка. 


Эффективность линеаризованной системы исследован путем сравнительного 
анализа в блочно-визуальной среде имитационного моделирования Simulink. 
Полученные положительные эффекты в задаче управления «двойного интегратора» 
позволить широко использовать результаты данной работы в инженерных приложениях. 


Ключевые слова: оптимальное управление, энергия управления, быстродействие, 
релейное управление, функция переключения, робастность, Simulink. 
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REDUCING THE ENERGY CONSUMPTION OF THE SPEED-OPTIMAL 
CONTROL SYSTEM 


Azerbaijan Technical University 


The article examines the influence of control vibration on energy consumption in a system 
that 1s optimal in terms of speed and the robust properties of a modified control system. The 
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solution of the classical problem of optimal control by the criterion of performance leads to 
two-position relay control. In this case, there are the following disadvantages: in a small 
neighborhood of the equilibrium position, vibration of the control signal occurs, which leads to 
excessive consumption of energy; the switching function is time dependent and has no 
feedback. Taking into account these shortcomings, various modifications of the problem 
formulations have been proposed. These statements differ by multi-criteria and. addition of 
functional restrictions. "Correction" of the original problem leads to the complication of the 
mathematical apparatus and the system-technical implementation of the system. One of the 
constructive approaches used in this work is the systems approach. In this case, the linearization 
of the optimal algorithm is carried out: the relay control is replaced by a linear-relay control, 
which in a small neighborhood of the equilibrium position allows you to connect feedback; the 
nonlinear switching function is approximated by a hyperplane, the equation of which includes 
state variables; in addition, the synthesis problem for high-order systems 1s simplified. 


The efficiency of the linearized system is investigated by comparative analysis in the 
block-visual simulation environment Simulink. The obtained positive effects in the control 
problem of a "double integrator" will make it possible to widely use the results of this work in 
engineering applications. 


Key words: optimal control, energy, speed, relay control, switching function, robustness, 
Simulink. 
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CƏLDİŞLƏMƏYƏ GÖRƏ OPTİMAL İDARƏETMƏ SİSTEMİNİN ENERJI 
SƏRFİNİN AZALDILMASI 


Azərbaycan Texniki Üniversiteti 


Məqalədə idarənin vibrasiyasının cəldişləməyə görə optimal idarəetmə sistemində enerji 
sərfinə təsiri tədqiq edilir. Klassik optimal idarəetmə məsələsinin cəldişləmə kriterisinə görö 
həlli ikimövqeli relt idarəsinə gətirir. Bu halda müəyyən çatışmamazlıqlar meydana çıxır: 
tarazlıq vəziyyətinin kiçik ətrafında idarə siqnalının vibrasiyası baş verir ki, bu da enerji sərfini 
artırır, qoşulma funksiuası əks əlaqəli olmayıb yalnız zamandan asılı olur. Göstərilən 
çatışmamazlıqları nəzərə alaraq hazırda məsələnin qoyuluşunun müxtəlif modifikasiyaları 
təklif edilmişdir. Bunlar çoxkriteriyalıq və müxtəlif funksional məhdudiyyətlərin əlavə 
olunması По fərqlənirlər. Bu yol ilo ilkin məsələnin “düzəlişi” riyazi aparatı və sistemin 
sistemotexniki realizasiyasını mürəkkəbləşdirir. Məqalədə konstruktiv üsullardan buri olan 
sistemli yanaşmadan istifadə edilmişdir. Bu halda optimal idarəetmə alqoritminin 
xəttiləşdirilməsi aparılır: rele idarəsi xətti-rele idarəsi ilə əvəz olunur ki, bu da tarazlıq 
vəziyyətinin kiçik ətrafında əks əlaqənin qoşulmasını təmin edir, qeyri-xətti qoşma funksiyası 
vəziyyət dəyişənlərindən asılı olan hipermüstovi ilo aproksimasiya olunur ki, bu halda yüksək 
tərtibli sistemlərin sintezi sadələşir. 
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Xottiosdirilmis sistemin effektivliyi Matlabin tərkibində olan vizual-bloklu imitasiya 
modelləşdirmə paketi Simulink-do müqayisəli analiz yolu ilə tədqiq edilmişdir. “İkiqat 
inteqrallay1cinin idarəetmə məsələsi”-ndə alınmış müsbət nəticələr işin nəticələrini mühəndis 
tətbiqlərində geniş istifadə etməyə yol açır. 


Açar sözlər: optimal idarəetmə, idarənin enerjisi, cəldişləmə, rele idarəsi, qoşma 
funksiyası, robastlıq, Simulink. 


1. Введение 


Задача оптимального управления по быстродействию в классической постановке 
приводит к увеличению затрат на реализацию оптимального управления. Основной 
причиной этой особенности является то что, оптимальное управление является 
релейным, принимающим предельные значения. При этом система имеет склонность к 
возникновению скользящего и автоколебательного в малой окрестности точки 
равновесия режимов. Наличие таких режимов приводит к перерасходу энергии и 
ограничивают срок службы исполнительных органов и других инерционных 
механизмов системы. 


Исходя из этого в [l] рассмотрена обобщенная задача с добавлением 
дополнительных ограничений в многокритериальной среде. 


В данной статье рассматривается альтернативная задача улучшения динамических 
и статических показателей системы путем модификации структурных особенностей 
классического решения. Задача быстродействия решена на примере управления 
«двойного интегратора». 


Объекты вида y —bu и задача двойного интегратора широко рассматривается в 


теории автоматического управления [2]. Двойной интегратор используется для 
моделирования электромеханических, электрических, гидравлических объектов. В 
частности, двигатель постоянного тока с независимым возбуждением |3,4,5 |, наведение 
самолета или ракеты на цель в простейшей задаче перехвата [6,7], мягкая посадка на 


Луну [8] и др. 


2. Исходная задача [9] 


В классической постанвке задача состоит из определения ограниченного 
управления — И < И < ЕИ , приводящего систему из заданного начального состояния 


х(0) = ху в начало координат Х(/,) = Ху = Оза кротчайшее время. Математическая 


постановка задачи быстродействия: 
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tr 
J = |а =, "HM | 
0 и<и 
x(t) = Ax(t) + Ви, 
-и<и<+и, 


х(0)=х, X(tp)=xX, =0. 


При отсутствии ограничения принимая и =оо, т.е. прикладывая к объекту 
бесконечную энергию, его можно было бы привести в точку равновесия мгновенно. 
Ясно, что при этом задача теряет физическую смысл. 


В этой задаче из-за наличия ограничения на управление определить решение из 


ОН 
необходимого условия оптимальности E — Ü не возможно. Оптимальное управление 
и 


определяется из условия тах{ Н}. Здесь функция Гамильтона 
Н =—1-- 4! (Ax + 4! (Ви. 


Оптимальное управление U(t) должно обеспечивать максимальное значение Н. 
Поскольку и входить только в третье слагаемое, то для достижения максимального 


значения функции Гамильтона должно выполняться условия: при es (t)>0 > из +и 


T - 
‚а при А (t) «0— и =—и. Таким образом оптимальное управление имеет 


следующий вид: 


и — usign(o), о- B'X(t). 


Полученное управление является управлением релейного типа. 


В общем случае сопряженные переменные определяются путем решения 
двухточечной задачи, состоящая из следующей канонической системы уравнений: 


1 
в= 1 = Ax Вили: n 
— — N 
öx 


x(0) + ху, xÜz) xr, 4(002 A9? 
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Непосредственное определение Х(/) и A(t) из решения уравнения (1) невозможно, 
поскольку начальное условие МО) = AQ является неизвестным. Поэтому при решении 
«параллельно» вычисляется и Ло, обеспечивающее попадание траектории x(t) в точку 


хт. В настоящее время разработаны различные алгоритмы вычисления Ао, в частности 
[10]. 

Для улучшения реализуемости и помехозащищенности системы, полученное B виде 
функции времени A(t) выражается в виде обратной связи A(X). Если для функции 
переключения с(х) невозможно получить аналитическое выражение или она 


получается слишком сложным, то прибегают к ее аппроксимации кусками более простых 
функций [11]. 


Не смотря на простоту оптимального управления и, = tu, = const оно не обладает 


стабилизирующим свойством. Для вырождения начало координат Ху = 0 в точку 
равновесия в ее малой окрестности управление должно принимать нулевое значение 
и, =O. Поскольку в составе оптимального управления нулевое составляющее 
отсутствует, то для удержания системы в малой окрестности начало координат 
управление совершает высокочастотную вибрацию с максимальной амплитудой. Такая 
структурная особенность приводит к перерасходу энергии. Поэтому при реализации 
релейного управления в малой окрестности начала координат необходимо включить 
нулевое управление или обратную связь. При этом можно использовать релейные 
регуляторы с статическими характеристиками, показанными на рис. |. 


Рис.1. Статические характеристики релейных регуляторов. 


а-двухпозиционное реле с зоной 
нечувствительностью; 


б- двухпозиционное реле с зоной линейностью 





(Saturation). 
a) 6) 


Пример 1. Рассмотрим задачу оптимального управления по быстродействию 
двойным интегратором. Уравнение объекта в координатах состояния имеет вид: 
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Требуется определить управление, приводящего объект из начального состояния 
Т s 
Ху = (Хо X29) в конечное x,-(0,0)' за наискорейшее время. 


Для данной задачи функция Гамильтона 


H=-1+),(t)x, +A, (tu. 


Оптимальное управление: 


и = мощ (t). 


Сопряженные переменные определяем из решения системы уравнений: 


ФАР | oH __ Т ә 29 
Fun A’ X(t) MO) = М =” 


Покоординатное представление этого уравнения: 


A (t) -0, A(t) = —Л. (1). 


При начальных условиях À (0) = A49, 42 (0) 229 решение: 


МО) = 4g, A(t) = Àog— МО. 


Тогда оптимальное управление до точности начальных условий: 


и = —sigh(A5g = Ма). 


В зависимости от Луи Аознак À (t) может изменится один раз. Не стараясь 
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находить начальные условия Лу, определим Л.(/) в виде обратной связи по 


состоянию О = Л„(х,Х,). 


Если определить уравнение фазовой траектории , TO при и = +! получим семейство 


2 2 
парабол х = —0.5x5 + С и X, = +0.5x5 + С. по которым изображающая точка двигается по 
1 2 1 1 2 2 


часовой стрелке (см. рис.3). Здесь Су, С, — зависящие от начальных условий постоянные. 


X. ^i 






Рис.3. Фазовый портрет оптимальной по 
G»0,u--] 





быстродействию системы 


Только два из парабол проходит через начало координат. 
Поэтому их можно использовать в качестве линий переключения 
управления. Уравнение линии, охватывающая обе части парабол можно записать в виде 


O =x, + 0.5хо | хо | = 0. 


Таким образом, оптимальное управление: 


-] if o»0, 
+1 Ё с<0. 


На рис.4 показана функциональная структурная схема системы с нелинейной 
функцией переключения с = 0.Ha рис.5,а показана переходная характеристика y(t), 0- 
сигнал управления U(t), в- фазовый портрет системы. 


ОБЪЕКТ 





x XY Graph 


Рис. 4 Функциональная схема оптимальной системы 
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a) , 6) 6) 


Рис. 5. Динамические характеристики системы 


Как видно из рисунка 5,6 в установивщем режиме сигнал управления совершает 
высокочастотное разрвное колебания (вибрация). 


3. Построение модифицированной системы 


В данной работе для устранения вибраций предложно структурные изменения 
оптимальной системы: нелинейная функция переключения O(X;, X») аппроксимируется 
прямой mei О; (Öz, Xə ) = x, +сх =0; sign- реле - заменяется на релейно-линейный 


закон управления (Saturation) . В общем виде предложенное управление: 


ks if is > So, 
u = sat(s) = 
usign(s) if is > $9. 


Ha рис.6 показаны статическая характеристика (a) и схема реализации (6) функции 
sat (-) на ракете Simulink. 





Gain Saturation 





а) 6) 
Рис.6. Статическая характеристика (а) и схема 


реализации (б) функции sat С.) 
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Здесь К, c - параметры настройки; $ = Ko, — входной сигнал; К =tg(a) = и / 54. 


В Simulinke задано К 21—» 5, = и. При этом на интервале линейности [s| < s) входной 
сигнал s проходит без изменения. Следовательно и = 5 = КО. 


Для данного примера с учетом $ = Ко, И и=1= 5 =1 выражение для функции 


sat(s) можно переписать в виде: 


и = sat(o,) << Ko, if о < бо, 


+1 if o, <о‹. 


Ha рис. 7 показана схема модифицированной системы, реализованной Ha Simulink 
c линейной функцией переключения o; = ху + сж =0,c=2. Коэффициент усиления 
K=20. 


Ha рис.8,а показана переходная характеристика y(t), б-сигнал управления u(t), 6- 
фазовый портрет системы. 


+-1 ОБЪЕКТ 


К 
> 5 | 1 у 


ain4 Saturation Transfer Fen? 


Isplay Integrator “Product 


Saturation‘ 


{| @ е(х1 ) 


Step 











2*u(1) 
В 
осоре! 


+-0.9 


XY Graph 


Puc.7. Схема модифицированной системы управления 
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a) 6) в) 


Рис.8. Динамические характеристики системы 


В установимщемся режиме сигнал управления принимает нулевое и(/)–0 значение и 
вибрации устраняются. Из puc.8, а и в видно, что время достижения начала координат 


(точка равновесия) составляет /, = 2c. 


4.Оценивание энергии управления и 
сравнительный анализ 


Энергию управления вычисляем по известному выражению 
T; 
2 
Q = |и 4, (2) 
0 


где Ís - время симуляции (время завершения процесса регулирования). 


В оптимальной системе время перехода t, = 2c. (см.рис.5 а,б). Энергия управления, 
вычисляемая по (2) при Т — 10c. составляет О=10 (см. рис.4, показание дисплея). В 
линеаризованной системе время перехода также t, = 2с. (см.рис.8, a, б). Однако при 


Г. =10с. > О =1.994 = 2 (см. рис.7, показание дисплея). 


Энергия управления оптимальной системы B 5 раз превышает энергии, 
затрачиваемой линеаризованной системой. Причина тому является вибрация 
управляющего сигнала в малой окрестности точки равновесия на интервале  f,«t«T; 
(рис.5,6). В системах стабилизации разница стремительно растет, поскольку в этих 
системах установившийся режим (settling time) может продолжаться довольно долгое 
время. 


Особый интерес представляет использование коэффициента усиления К. В 
определенных условиях увеличение К приводит к улучшению робастных свойств 
системы. 
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5. Исследование робастности 


Если расширить ограничение на управления, предложенный алгоритм можно 


использовать для управления объектов с параметрическими неопределенностями 
Покажем эту особенность на конкретном примере. 


Пример 2. Рассмотрим гармонический осциллятор с передаточной функцией 
W = b/(sf” +a). 


Параметр а изменяется на широком интервале 1<a<10 


На рис.9 показана схема управления с эталонной моделью g=sin(2t) (уставка) 


Ограничение на управление не учитывается. Это достигнуто путем устранения блока 
«Saturation». 





g=sin(4t) 


5 ок | B p 
Sine Wave т | 


In Н In ator 1 
Gain? E iar 











и 


Баһвжабоп1 


Рис.9. Схема системы управления при параметрической 


неопределенности 





Rustamov R. Q: REDUCING THE ENERGY CONSUMPTION OF... 72 





Рис.10. Динамические характеристики системы управления при 


параметрической неопределенности 


Как видно из рис.10,а при К=20 после t = 1.8 с. выход системы (у (Г), yo (t), уз (®)} 
довольно точно отслеживает эталонную траекторию g —sin( 27). Ощутимый разброс 


траекторий в установившемся режиме не наблюдается. Сигнал управления изменяется 


на небольшом интервале S<u,<5, 1=1,2,3. Полученный результат подтверждает 


наличие свойств робастности. Увеличивая коэффициента усиления К ошибку 
регулирования можно уменьшить до сколь угодно малой величины. 


6. Заключение 


Определение управления по быстродействию путем решения задачи аналитического 
конструирования оптимальных регуляторов является сложной математической задачей. 
Кроме того, полученное решение не является полным, поскольку начальное условие 
сопряженных переменных является неизвестным. В статье предложен структурный 
подход, состоящий из линеаризации оптимального релейного управления и функции 
переключения. При линеаризации использован принцип обратной связи. Основной 
задачей являлся не сохранение полной оптимальности, а уменьшение энергии 
управления путем устранения вибрации и наделения системы свойством робастности. 


Эффективность предложенного алгоритма модификации структуры проверена путем 
сравнительного анализа на Simulink. Полученные положительные результаты в задаче 
управления «двойного интегратора» позволить использовать результаты данной работы 
в инженерных приложениях. 
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